Lectures on the integrability of the 6-vertex 

model. 

Nicolai Reshetikhin * 

Contents 

1 Introduction Q 

2 Classical integrable spin chains [H 

2.1 Classical r-matrices and the construction of classical integrable 
spin chains 

2.2 Classical L-operators related to s?2 

2.3 Real forms 

2.3.1 

2.3.2 

2.4 Integrable classical local SC^-spin chains 



3 Quantization of local integrable spin chains |k] 

3.1 Quantum integrable spin chains 

3.2 The Yang-Baxter equation and the quantization 

3.3 Quantum Lax operators and representation theory 

3.3.1 Quantum LSL2 

3.4 Irreducible representations 

3.4.1 

3.4.2 

3.4.3 

3.4.4 

3.4.5 

3.4.6 Real forms 

3.5 The fusion of R- matrices and the degeneration of tensor products 
of irreducibles 

3.6 Higher transfer-matrices 

3.7 Local integrable quantum spin Hamiltonians 

3.7.1 Homogeneous ST/ (2) spin chains 



'Department of Mathematics, University of California at Berkeley, Berkeley, CA 94720- 
3840. and KDV Institute for Mathematics, Universiteit van Amsterdam, Plantage Muider- 
gracht 24, 1018 TV, Amsterdam, The Netherlands. E-mail: reshetik@math.berkeley.edu 



1 



3.7.2 Inhomogeneous ST/ (2) spin chains \2£ 

The spectrum of transfer-matrices 

4.1 Diagonalizability of transfer-matrices 

4.2 Bethe ansatz for sl% 

4.3 The completeness of Bethe vectors 

4.3.1 

4.3.2 

4.3.3 



5 The thermodynamical limit 

6 The 6-vertex model 

6.1 The 6-vertex configurations and boundary conditions 

6.2 Boundary conditions 

6.3 The partition function and local correlation functions 

6.4 Transfer-matrices 

6.5 The commutativity of transfer-matrices and positivity of weights 

6.5.1 Commutativity of transfer-matrices 

6.5.2 The parametrization 

6.5.3 Topological nature of the partition function of the 6-vertex 
model 

6.5.4 Inhomogeneous models with commuting transfer-matrices 

6.5.5 

6.5.6 

7 The 6-vertex model on a torus in the thermodynamic limit 

7.1 The thermodynamic limit of the 6-vertex model for the periodic 
boundary conditions 

7.2 The large N limit of the eigenvalues of the transfer-matrix .... 

7.3 Modularity 

7.3.1 

7.3.2 



8 The 6-vertex model at the free fermionic point 

8.1 Homogeneous case 

8.1.1 

8.1.2 

8.2 Horizontally inhomogeneous case 

8.3 Vertically non-homogeneous case 

8.3.1 

8.3.2 Largest eigenvalue 



2 



9 The free energy of the 6-vertex model 

9.1 The phase diagram for A > 1 

9.1.1 

9.1.2 

9.1.3 

9.2 The phase diagram |A| < 1 

9.3 The phase diagram A < — 1 

9.3.1 The phase diagram 

9.3.2 The antiferromagnetic region 

10 Some asymptotics of the free energy 

10.1 The scaling near the boundary of the D-region 

10.2 The scaling in the tentacle 

10.3 The 5- vertex limit 

10.4 The asymptotic of the free energy near the tricritical point in the 
5- vertex model 

10.5 The limit A -> -1" 

10.6 The convexity of the free energy 

11 The Legendre Transform of the Free Energy 

11.1 

11.2 

12 The limit shape phenomenon 

12.1 The Height Function for the 6-vertex model 

12.2 Stabilizing Fixed Boundary Conditions 

12.3 The Variational Principle 

12.3.1 

12.3.2 The variational principle 

12.3.3 

12.4 Limit shapes for inhomogeneous models 

12.5 Higher spin 6-vertex model 

13 Semiclassical limits 

13.1 The semiclassical limits in Bethe states 

13.2 The semiclassical limits in the higher spin 6-vertex model . . . . 

13.3 The large N limit in spin- 1/2 spin chain as a semiclassical limit . 

13.3.1 

13.3.2 

13.3.3 

14 The free fermionic point and dimer models 



3 



A Symplectic and Poisson manifolds 

A.l 

A.2 

A.3 



B Classical integrable systems and their quantization |75 

B. l Integrable systems in Hamiltonian mechanics ]75 

C Poisson Lie groups 

C. l Poisson Lie groups 

C.2 Factorizable Poisson Lie groups and classical r-matrices 

C.3 Basic example G = LSL2 

C.4 Symplectic leaves 



D Quantization 

D.l Quantization 

D.l.l Quantized algebra of observables 
D.2 Examples of family deformations .... 

D.2.1 

D.2.2 

D.3 Quantization of integrable systems . . . 



1 Introduction 

The goal of these notes is to outline the relation between solvable models in 
statistical mechanics, and classical and quantum integrable spin chains. The 
main examples are the 6- vertex model in statistical mechanics, and spin chains 
related to the loop algebra Ls/ 2 . 

The 6-vertex model emerged as a version of the Pauling's ice model, more 
generally, as a two-dimensional model of ferroelectricity. The free energy per 
site was computed exactly in the thermodynamical limit by E. Lieb |23j . The 
free energy as a function of electric fields was computed by Sutherland and Yang 
|37] . For details and more references on earlier works on the 6-vertex model see 
[24] . The structure of the free energy as a function of electric fields was also 
studied in [27][B]. 

Baxter descovered that Boltzman weights of the 6-vertex model can be ar- 
ranged into a matrix which staisfies what is now known as the Yang-Baxter 
equation. For more references on the 6-vertex model and on the consequences 
of the Yang-Baxter equation for the weights of the 6-vertex model see [3] . 

The 6-vertex model and similar 'integrable' models in statistical mechanics 
became the subject of renewed research activity after the discovery of Sklyanin 
of the relation between the Yang-Baxter relation in the 6-vertex model and the 
quantization of classical integrable systems [M] . It lead to the discovery of many 
'hidden' algebraic structures of the 6-vertex model and to the construction of 
quantizations of a number of important classical field theories [13] [14] . For more 
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references see, for example [H] . Further development of this subject resulted in 
the development of quantum groups [llj and further understanding of algebraic 
nature of integrability. 

Classical spin chains related to the Lie group SL2 is an important family 
of classical integrable systems related to the 6-vertex model. The continuum 
version of this model is known as Landau-Lifshitz model. One of its quantum 
counterparts, the Heisenberg spin chain has a long history. Eigenvectors of its 
Hamiltonian of the Heisenberg model were constructed by Bethe in 1931 us- 
ing the substitution which is known now as the Bethe ansatz. He expressed 
the eigenvalues of the Hamiltonian in terms of solutions to a system of alge- 
braic equations known now as Bethe equations. An algebraic version of this 
substitution was found in |14) . 

These lectures consist of three major parts. First part is a survey of some 
basic facts about classical integrable spin chains. The second part is a survey of 
the corresponding quantum spin chains. The third part is focus on the 6-vertex 
model and on the limit shape phenomenon. The Appendix has a number of 
random useful facts. 

The author is happy to thank A. Okounkov, K. Palamarchuk, E. Sklyanin, 
and F. Smirnov for discussions, B. Sturmfels for an important remark about the 
solutions to Bethe equations, and the organizers of the school for the opportunity 
to give these lectures. 

This work was supported by the Danish National Research Foundation through 
the Niels Bohr initiative, the author is grateful to the Aarhus University for the 
hospitality. It was also supported by the NSF grant DMS-0601912. 

2 Classical integrable spin chains 

The notion of a Poisson Lie group developed from the study of integrable systems 
and their relation to solvable models in statistical mechanics. 

It is a geometrical structure behind Poisson structures on Lax operators, 
which emerged in the analysis of Hamiltonian structures in integrable partial 
differential equations. First examples of these structures are related to taking 
semiclassical limit of Baxter's i?-matrix for the 8-vertex model. The notion of 
Poisson Lie groups and Lie bialgebras places these examples into the context 
of Lie theory and provides a natural versions of such systems related to other 
simple Lie algebras. 

2.1 Classical r-matrices and the construction of classical 
integrable spin chains 

Here we will recall the construction of classical integrable systems based on 
classical r-matrices. 

A classical r-matrix with (an additive) spectral parameter is a holomorphic 
function on C with values in End(V)® 2 which satisfies the classical Yang-Baxter 
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equation: 



[ri 2 {u),r 13 (u + v)} + [ri 2 (u),r 23 {v)} + [r 13 (u + v), r 2 3(v)] = (1) 

where rij(u) act in U® 3 , such that Ti 2 {u) = r(u) g3 1, r 23 (u) — 1 <g) r(u), etc. 

Let L(u) be a holomorphic function on C of certain type (for example a 
polynomial). Matrix elements of coefficients of this function satisfy quadratic 
r- matrix Poisson brackets if Their generating function L(u) has the Poisson 
brackets 

{Li(u),L 2 (v)} = [r(u),Li(u)L 2 (v)] (2) 

where Li{u) = L{u) (g> l,Z<2(w) — 1 <S) L(u). The expression on the left is the 
collection of Poisson brackets {Lij(u), Lki(v)}. 
Consider the product 

T(u) = L( N \u-a N )...L { i\u-ai) 

Matrix elements of T(u) are functions on P^ x • • • x P x . The factor Lfo is the 
function on The Poisson structure on Pj is as above. 

The r- matrix Poisson brackets for L^(u) imply similar Poisson brackets for 
Tij(u): 

{T 1 («),r 2 (w)} = [r(u-v),Ti(u)T 2 (v)] 
Taking the trace in this formula we se that 

{t(«),t(«)}=0 

Fix symplectic leaves Si C Pi for each i.The restriction of the generating 
function t(u) to the product of symplectic leaves S = Sn x ■ ■ • x Si C Pn x 
■ ■ • x Pi gives the generating function for commuting functions on this symplectic 
manifold. 

Under the right circumstances the generating function t(u) will produce the 
necessary number of independent functions to produce a completely integrable 
system, i.e. dim{S)/2. 

One of the main conceptual questions in this construction is: Why classical 
r-matrices exist? Indeed, if n = dim(V), the equation (fTJ) is a system of n 6 
functional equations for n 4 functions. The construction of the Drinfeld double 
of a Lie bialgebra provides an answer to this question [TT] . 

2.2 Classical L-operators related to s/ 2 

In this section we will focus on classical L-operators for the classical r-matrix 
corresponding to the standard Lie bialgebra structure on sl 2 . 

Such L-operators describe finite dimensional Poisson submanifolds in the 
infinite dimensional Poisson Lie group LSL 2 . For some basic facts and refer- 
ences see Appendix [X] Up to a scalar multiple they are polynomials in the 



G 



spectral variable z. One of such simplest Poisson submanifolds correspond to 
polynomials of first degree of the following form: 

where z = exp(u). The r-matrix Poisson brackets on L(u) with the r-matrix 
given by (f9"8j) induce a Poisson algebra structure on the algebra of polynomials 
in a,b, .... This Poisson algebra can be specialized further (by quoticnting with 
respect to corresponding Poisson ideal). As a result we arrive to the following 
L-operator: 

which satisfies the r-matrix Poisson brackets @ with the following brackets on 
k,e,f: 

{k,e} = eke, {k,f} = -ekf (5) 
{e,/} = 2e(fc 2 -fc- 2 ) ; (6) 

The function 

c = ef + k 2 + k- 2 , (7) 

Poisson commute with all other elements of this Poisson algebra, i.e. it is the 
Casimir function. It is easy to show that this is the only Casimir function on 
this Poisson manifold. 
It is easy to check that 

L(-tt)' = -D z olL{u)alD- 1 (8) 

On the level surface of c parameterized as c = t 2 + t~ 2 this matrix satisfies 
the extra identity 

L(u)6(L(~u)Y = (zt - x-H- 1 )^- 1 - t _1 z)I 

where I is the identity matrix and 

9(e) = /, 9(f) = e, 9(k) = k 

is the anti-Poisson involution: {9(a)9(b)} = — 9({a, b}), 9 2 = id. 
Its easy to find the determinant of L(u): 

det(L(u)) = (zt - z^ x r x )(r l z - tz^ 1 ) 

When z — t, t^ 1 the matrix L degenerates to one dimensional projectors. 

T=TZ \ o /+-!„ +1,-1 



L(t) = ( -TJe J <g, (t~i ej tk- 1 - t~ l k) (9) 

/ t~ 1 k-tk~ 1 \ 

L(t~ 1 )=l te J (g> (te.t^k- 1 -tk) (10) 
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2.3 Real forms 



2.3.1 

Here we will describe the Poisson manifold SU2 ■ Let Si , S2 , S3 be the coordi- 
nates SU2 corresponding to the usual orthonormal basis in su2- The Poisson 
brackets between these coordinate functions are 

{Si, S2} = {S2, S3} = 2Si, {S3, Si} = 25*2, 

These coordinates are also known as classical spin coordinates 
this Poisson algebra is generated by C = S 2 + Sf + Sf. 

It is convenient to introduce S + = (Si + iS 2 )/2, S~ = (Si 
Si are real, S + = S~ . 

The Poisson brackets between these coordinates are: 

{S + ,S-} =-iS s , {S 3 ,S ± }=t2iS ± 

Level surfaces of C are spheres. On the level surface with C = I 2 with two 
point S3 = ±1 being removed we have the following Darboux coordinates: 

S+ = yj pi - p* , S~ = e-'ty/pl-p 2 , S 3 = l-2p 
where < p < I and < <f> < 2tt and {p, <fr} = 1. 

2.3.2 

The Poisson algebra 

{k,e} = eke, {k,f} = -2ekf, 
{e,f} = 2e(k 2 -k- 2 ) 

has the two real forms which are important for spin chains with compact phase 
spaces: 

• The real form with e = 1, |fc| = 1, e = /. 

• And the real form where e = i, e = / , and k = k. 

In the first case the level surface of (|7|) with c = 2 cos 2R without two points 
e = / = has the following Darboux coordinates: 

e = 2e** Vsin(p) sin(2i? - p), f = 2e -i Vsin(p) sm ( 2 # - p), k = e^ R ~^ 

with {p, 4>} = 1, and < p < I and < 4> < 2n. 

Similarly, in the second case the level surface c = 2 cosh R without two points 
e = / = have Darboux coordinates 

e = 2e^ Vsinh(p) sinh(2i? - p), f = 2e~^ ^&mh(p) sinh(2i? - p), k = e R ~ p 

where 4> G R, < p < R, and {p, 0} = 1. 

We will denote these level surfaces by S^ R '. In the compact case S^ is dif- 
feomorphic to a sphere. It can be realized as the unit sphere with the symplectic 
form dependent on R. 



. The center of 
~iS 2 )/2. Since 
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2.4 Integrable classical local S^-spin chains 

Because the S^ R ' is a symplectic leaf of the Poisson Lie group LGL 2 , the Carte- 
sian product 

m R u -,Rn = S (Ri) x ... x s(-Rjv) 

is a Poisson submanifold in LGL 2 , which means matrix elements of the mon- 
odromy matrix 

T(z)^L^\za 1 )L (R ^(za 2 )...L (RN Hza N ) (11) 

satisfy the r-matrix Poisson brackets. 

In order to obtain local Hamiltonians in the homogeneous classical spin chain 
ai = ■ ■ ■ = O/v = 1, i?i = • • • = Rjy = R one can use the degenerations ©, ([TO)) . 
Combining these formulae ([9]), (ITO| . and (|8j we obtain the following identities: 

L(t)=a®P\ L(t- 1 ) = -o v DP®a t D- 1 a y 

where column vector a and row vector j3 are given in (|9]). These identities 
imply 

N N 

tr(T(t)) = H(p n ,a n+1 ), triTit- 1 )) = (-l) N 

71—1 71 — 1 

Here we assume the periodicity a_/v+i = ai and (3n+i = Al- Now notice that 

One the other hand this trace can be computed explicitly: 

tr(L n {t)L n+ i(t)) = e n f n+ i + f n e n +i + (tk„ - t _1 fc^ 1 )(ifc„ + i - 

(tk- 1 - r^Kt^j - r x fc n+1 ) (12) 

This gives the first local Hamiltonian 

N 

H = log(tr(T(t))tr(T(r 1 ))) = £ log(e n / n+1 + /„e n+1 + 

n=l 

(tk n - t^k-^itkn+i - r 1 ^) + (tk- 1 - rXX^ji - rX+O) (is) 

Other local spin Hamiltonians can be chosen as logarithmic derivatives of 
tr(T(z)) at when z = t ±x , for details see P3] and references therein. 

When N is even and inhomogeneities are alternating ai — a,a 2 — fl , 03 = 
a, . . . , cl 7v = a" 1 there is a similar construction of local Hamiltinians also based 
on degenerations (JTOJ. Again, all logarithmic derivatives of T(z) at points 
z = at^ 1 ,a _1 i ±:L are local spin Hamiltonians. 

In the continuum limit the Hamiltian dynamics generated by these Hamil- 
tians converges to the Landau-Lifshitz equation, see |14j and references therein. 
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3 Quantization of local integrable spin chains 



3.1 Quantum integrable spin chains 

In the appendix ?? there is a short discussion of integrable quantization of 
classical integrable system. 

A quantization of a local classical integrable spin chin is an integrable quanti- 
zation of a classical local integrable spin chain such that quantized Hamiltonians 
remain local. That is the collection of the following data: 

• A choice of the quantization of the algebra of observables of the classical 
system (in a sense of the section Appendix ??), i.e a family of associative 
algebras with a ^-involution which deform the classical Poisson algebra of 
observables. 

• A choice of a maximal commutative subalgebra in the algebra which is a 
quantization of Poisson commuting algebra of classical integrals. 

• In addition, locality of the quantization means that the quantized algebra 
of observables is the tensor product of local algebras (one for each site 
of our one-dimensional lattice): = <8>^ r =1 B^ n \ and that the quantum 
Hamiltonian has the same local structure as the classical Hamiltonian 
(??): _ 

n 

where H n = 1®- • -(g)ifW(g)- • -®1 and H W e B < ^ ) ®B < ^ +1) (E)---(E)B^ +k) . 

• A ^-representation of the algebra of observables (the space of pure states 
of the system). 

3.2 The Yang-Baxter equation and the quantization 

Here we will describe the approach to the quantization of classical spin chains 
with r-matrix Poisson bracket for polynomial Lax matrices based on construc- 
tion of corresponding quantum i?-matrices and quantum Lax matrices. 

In modern language the construction of quantum i?-matrix means the con- 
struction of the corresponding quantum group, and the construction of the quan- 
tum L-matrix means the construction of the corresponding representation of the 
quantum group. 

Suppose we have a classical integrable system with commuting integrals 
obtained as coefficients of the generating function t(u) described in section I2TT1 
The i?-matrix quantization means the following: 

• Find a family R(u, h) of invertible linear operators acting in V <E> V such 
that for each h they satisfy the quantum Yang-Baxter equation 

i?i 2 (u, h)Ri 3 (u + v, h)R 2 3{v, h) = i?2s(«, h)Ri 3 (u + v, h)Ri 2 (u, h) 
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and when h — > 

R(u,h) = 1 + hr{u) + 0(h 2 ) 
where r{u) is the classical r-matrix. 

• Let the classical Lax matrix L(u) be a matrix valued function of u of cer- 
tain type (for example a polynomial of fixed degree) , with matrix elements 
generating a Poisson algebra with Poisson brackets @. For given R(u, h) 
define the quantization of this Poisson algebra as the associative algebra 
generated by matrix elements of the matrix £(u) of the same type as L{u) 
(for example a polynomial of the same degree) with defining relations 

R(u, h)L(u + v) <®L{v) = (l® L(v)){L(u + v) ® l)R(u, h), (14) 

Denote such algebra by B^. 

• Consider the generating function 

T{u) — L\{u — wi) ® . . .Ln(u — wn) (15) 

acting in End{V)®B^ ® ■ ■ -QB^. It is easy to see that the commutation 
relations (|14p imply 

R(u, h)T(u + v)® T(v) = (1 (g) T(v))(T(u + v) ® l)J?(u, /i), (16) 

The invertibility of R(u,h) together with the relations (|16|) imply that 
t(u) — trv(T(u)) £ B^ ® • •• ® is a generating function for a com- 
mutative subalgebra in i?^ 1 ' ® • •• ® B^: 

[t(u),t(v)]=0 

Under the right circumstances this commutative subalgebra is maximal 
and defines an integrable quantization of the corresponding classical inte- 
grable spin chain. 

The i?-matrix was found by Baxter. Sklyanin discovered that when h — > 
the classical i?-matrix defines a Possion structure on LGL2 defined by the 
formula (|100[) and that it implies the Poisson commutativity of traces. 

There is an algebraic way to derive the Baxter's R-matrix from the universal 
i?-matrix for Uqijglz)- It is outlined in the appendix. 

3.3 Quantum Lax operators and representation theory 
3.3.1 Quantum LSL 2 

Here is the formal definition of C q (SL2) in terms of generators and relations. 
Let q be a nonzero complex number. The algebra C q (SL2) is a complex algebra 
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generated by elements the where i,j = 1,2 and k e Z. Consider the 

-(fc) 



matrix 7{z) which is the generating function for the elements 



°° / rp(Q 

k=i V 



(0) T (o) \ 

1 22 J 



The determining relations in C q {SL 2 ) can be written as the following matrix 
identity with entries in C q {SL 2 ): 

R(z)7(zw)®7(w) = (l®7(w))(7(zw)®l)R(z), (18) 
7{qz) n 7(z) 22 - 7(qz) 12 7(z) 21 = 1, 

where the tensor product is the tensor product of matrices. Matrix elements in 
this formula are multiplied as elements of C q (GL 2 ) in the order in which they 
appear. 

The matrix R(z) acts in C 2 ®C 2 and has the following structure in the tensor 
product basis: 

/l 0\ 

f(z) z~ l g{z) 
zg(z) f(z) 
\0 1/ 

where 



i?(z) 



(19) 



ft \ z-z 1 (q-q x ) 

/(*) = rr-ri' 9 Z = rr— T 

zq — z L q L zq — z L q 1 

It satisfies the Yang-Baxter equation. 



Remark 1. There is important function s{z): 
s(z)=q- 1 / 2 



I 2 2 4\2 



where 



(*W)oo(*V;5 4 )c 

OO 

^;p)oc=n( 1 -^")- 



i/ie matrix 



R(z) = s(z)R(z) 

satisfies what is known in physics unitarity and the crossing symmetry: 
^(z)^(z- 1 )* = 1, ^(z- 1 )' 2 = C^zq- 1 )^ 

where t is the transposition with respect to the standard scalar product in C 2 ®, 
t 2 is the transposition with respect to second factor in the tensor product, and 
C 2 = 1 ® C where 
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The Hopf algebra structure on C q {SLi2) is determined by the following action 
of the antipode on the generators: 

A(T tf (z))= 7 ik {z)®7 kj {z). (20) 

fc=l,2 

The right side here, as well as in the second relation in (fl8l) , is understood as a 
product of Laurent power series. The antipode is determined by the relation 

7{z){S® id7(z)) = 1. 

Let d be a nonzero complex number. The identity (|T8j) implies that the 
power series 

T 1 {z;d)=d7 11 {z)+d- 1 7 22 {z) (21) 

generates a commutative subalgebra in C q (SL2). 

Chose a linear basis (for example ordered monomials in ) . The relations 
between generators will give the multiplication rule for monomials which will 
depend on q. This multiplication turns into the commutative multiplication of 
coordinate functions when q = 1. The commutator of two monomials, divided 
by q — 1, at q = 1 becomes the Poisson bracket. It is easy to check that this 
Poisson brackets is exactly the one defined by the classical r-matrix (j98| . 

Remark 2. Let D be a diagonal matrix. It is easy to see that [D <g> 1 + 1 ® 
D, R(x)] =0. It is easy to show that if 

R(z)7 1 (zw)7 2 (w) =7 2 (w)7x(zw)R(z) 

then 

R(z) = (z D ® l)R(z)(z- D ® 1), T(z) = z D 7(z)z- D 
satisfy the same relation 

R{z)7x{zw)72{w) =7 2 (w)7i(zw)R(z) 

In particular R(z) satisfies the Yang- Baxter equation. Choosing D — diag(— 1/2,1/2) 
gives the R-matrix hi 9(1 but with no factors z ±x off-diagonal. This symmetric 
version of the R-matrix is the matrix of Boltzmann weights in the 6-vertex model. 

Remark 3. If A is an invertible diagonal matrix such that (A ® A)R(z) = 
R(z)(A ® A) and7(z) is as above then 

7 A {z) = A7(z)A- 1 

also satisfies the relations H8\ ). 
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3.4 Irreducible representations 
3.4.1 

It is easy to check that the following matrix satisfies the i?-matrix commutation 
relations (|T51) 

z(z) = ( I***-*-; 1 *- 1 ?-* S 1<ri l * ) (22) 

\ zq 2 e zk —z kq 2 / 

if e, f,k commute as 

ke = qek, kf — q~ x fk, 
ef-fe = (q-q-^e-k- 2 ). 

Denote this algebra C q . 
The element 

c = fe + ^q + k^q- 1 (23) 

generates the center of this algebra. 

It is clear that this algebra quantizes the Poisson algebra (O©. Indeed, 
the algebra C\ is the commutative algebra generated by e,/, fc ±1 . Consider 
the monomial basis e n k m f in G q . Fix the isomorphism between C q and C\ 
identifying these bases. The associative multiplication in C q is given in this 
basis the function of q: 

k 

it is clear that when q = 1 this multiplication is the usual multiplication in the 
commutative algebra generated by e, /, fc, The skew symmetric part of the 
derivative of m(q) at q = 1 is the Poisson structure. 

The algebra C q is a Hopf algebra with the comultiplication acting on gener- 
ators as 

Ak = k ® k, Ae = e ® k + k' 1 ® e, A/ = / ® + A; -1 ® / 

The algebra C g is closely related to the quantized universal enveloping alge- 
bra for sl2- Indeed, elements E = ekj (q — q^ 1 ), F = k^ 1 f /(q — q^ 1 ), K = k 2 
are generators for U q (sl2)- 

K — K^ 1 

KE = q 2 EK, KF = q~ 2 FK, EF - FE = — 

q-q 

with 

AK = K ® K, AE = E <E> K + 1<E) E,AF = F ®1 + K^ 1 <Z> F 
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3.4.2 

Assume that q is generic. Denote by the irreducible m + 1-dimensional 

representation of Cq, and by Vq the highest weight vector in this representa- 
tion: 

bj m ' = q^VQ m \ ev^ = 

The weight basis in this representation can be obtained by acting / on the 
highest weight vector. Properly normalized the action of C q on the weight basis 
is: 

kv^ = q^- n V^\ M m > = (q m - n - q- m+n )v%l\, ev™ = (q n - q~ n )v^\. 

The Casimir element c acts on V(m) by the multiplication on q m+1 +q~ m ~ 1 . 
Because the algebra C q is a Hopf algebra, it acts naturally on the tensor 
product of representations. 

3.4.3 

Denote the matrix (|22[) evaluated in the irreducible representation l/( m ) by 
£( m )(z). It is easy to check that it satisfies the following identities: 

£ (m)( z -i)t = -D z CL {m) {zq)C- l D- 1 

L^izfL^iz- 1 ) = {zq^ - z- 1 q-^)( Z - 1 q' Ii f 1 - zq~^)I 

where D z is a diagonal matrix. 

Here t is the transposition with respect to the standard scalar product in 
C 2 and T is the transposition t combined with the transposition in T/( m ) with 
respect to the scalar product , tij^J 1 ) = 5 n ,ri- 

3.4.4 

The matrix (|22|) defines a family of 2-dimensional representations of C q (SL2). 
If a is a non-zero complex number such representation is 

7(z) ^ g m (z)^ m \za) 

where the factor g m {z) is important only if we want to satisfy the second re- 
lations which play the role of quantum counter-parts of the the unimodularity 
(i.e. det = 1) of TO). 

M __ rW;g 4 )oo(^ 2 ;g 4 )co 
9m[z> *(*-V 2 ;<? 4 )oo(*9 3 *W)oo 
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3.4.5 



The comultiplication defines the tensor product of irreducible representations 
described above 

N 

T(z) II ffm*(«/oi)T (,ni '- ,ffw) («|oi. • • • > «n) (24) 
i=i 

where 

T^-> m »\z\ ai , ...,a N ) = L x {z/ai) . . . JC N (z/a N ), (25) 

where we have taken the matrix product of the matrices £(2) as 2 x 2 matrices. 
The n-th factor in fl55j} acts on the rc-th factor of <g) • • • ® 0™"). 

Notice that this is also a tensor product of representations of C g . 

3.4.6 Real forms 

As in the classical case there are two real form of the algebra C q which are 
important for finite dimensional spin chains. 

Recall that a *-involution of a complex associative algebra is an ani- involution 
of the algebra, i.e. (ab)* = b*a* which is complex anti-linear: (Aa)* = Aa*. A 
real form of a complex A corresponding to this involutions is real algebra which 
is the real subspace in A spanned by the ^-invariant elements. 

When \q\ — 1, we will write q = exp(i7), the relevant real form of C q , is 
characterized by the ^-involution which acts on generators as 

e* = /, k* = k- 1 

When q is real positive we will write q = exp(jj). In this case the relevant 
real form is characterized by the * involution acting on generators as: 

e = ,/, / = e, k =k 

As 7 — > or r\ — > these real forms become real forms of the corresponding 
Poisson algebras described in section 12.31 with e = i and e = 1 respectively. 

3.5 The fusion of i?-matrices and the degeneration of ten- 
sor products of irreducibles 

This section is the analog of the construction of quantum L-operators by taking 
the tensor product of 2-dimensional representations. 

Consider the product of i?-matrices acting in the tensor product of n + m 
copies of C 2 : 

Rl>2'...m',12...n{z)= Rl'l{z) R V2 (zq) ... i?i'„(zq" _1 ) 

Rvi(zq) R 2 >2(zq 2 ) ... R 2 >n(zq n ) 

Rmnizq™- 1 ) Rm>2(zq m ) ... R m > n {zq n + m - 2 ) 
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The operator R(z) satisfies the identities 

R{z)R t (z- 1 ) = (zq - z^q-^iz^q - zq' 1 ) 
PR(z)P = R t (z) 

where t is the transpiration operation (with respect to the tensor product of the 
standard basis in C 2 ), and 

det(R(z)) = (zq - z- 1 <r 1 ) 3 0~ 1 <? - m~ X ) 

From here we conclude that the matrix PR(z) degenerates at z — q and z = q^ 1 
to the matrices of rank 3 and 1 respectively. 
Define 

Pu n = Ria(«) R 23 (g 2 ) ... Rmiq 71 - 1 ) 

R 23 (<z) •■• R2n(q n - 2 ) (26) 

Rn-lnW 

where R(z) = PR(z) and P is the permutation matrix, P(x ® y) = y ® x. 

Consider (C 2 )®" as a representation of C q . Because all finite dimensional 
representations of this algebra are completely reducible, it decomposes into the 
direct sum of irreducible components. The irreducible representation y( n ) ap- 
pears in this decomposition with multiplicity 1. One can show that the operator 
(|26|) is the orthogonal projector to V^ n \ The proof can be found in |22) . 

Also, it is not difficult to show that 



12. ..n 



(zq-^ - z~ l q^) . . . (zq^ - 2 _1 ff _ V )^L.n](«) ( 27 ) 



where the linear operator R 1,n (z) acts in C 2 ®^™* 1 and the second factor appears 
as the g-symmetrized part of the tensor product C 2 ® . Moreover, it is easy to 
show that this operators is conjugate by a diagonal matrix to &(z): 

R^ l ' n \z) : 



zkqz — z 1 k 1 q 2 z x q 

zqze zk~ l q? — z~ 1 kq~~; 

where e, /, k act in the n + 1 dimensional irreducible representation as it is 
described in section 13.41 In this realization of the irreducible representation 
weight vectors appear as v^ 1 ' = P±,_ n ei ® e± <g> e-i <g> e-i where we have n — k 
copies of ei and k copies of e2 in this tensor product. 
Similarly 

PU...m'Pl2...n R ^^ZZZJ ■•• Rl ' n ( Z €T + l 

Ryi{zq~ n * ) ... R2'n(zq" ™ ) 



p+ p 

r V2'...m< -M 



12. ..n 
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(zq- 
(zq- 



2 — Z 1 q 2 J 



(zq 2 — z 1 q 



z 1 q 



(zq—^ 



z 1 q 



#( m '")(z) 

Here we assume that m < n. The matrix elements of R^ m ' n \z) are Laurent 
polynomials of the form z~ m P(z ) where P(t) is a polynomial of degree m. 

The matrix R^ m ) also can be expressed in terms e, /, k . 

The matrices RS k < l > satisfy the the Yang-Baxter equation 

R^\z)R[ l 3 n \zw)R^ n \w) = R%> n \w)R^ l \z W )R[ l 2 m) (z) 
In addition to this they satisfy identities 



R^ m \z)R i - m ^{z- l f 



2i(z)s ml (z x ) 



and 



R^\z)^ = (-l) l (D z CW ® ^R^izq^D- 1 ^ 1 ® 1) 
where s m i(z) is a Laurent polynomial in z which is easy to compute, C^ m ^ = 
Py2 n ® CPj^ m an d we assume that I < m. 

3.6 Higher transfer-matrices 

Define C^SI^) -valued matrices 

where P^ 2 n is defined above and T is the matrix (fl7|) . 
They satisfy the relations 

R (Lm) (z) 12 l[ l \zw)J { 2 m) (w) = J i 2 m) (w)J i 1 l) (zw)R^ m \z) 12 

For non-zero d define the following elements of C q {SL 2 ) 

r e (z) = (id ®tr v(e) ){7^{z)d {l) ) . 

where = diag(d l ,d l ~ 2 , . . .,d~ l ) and d ^ 0. 

The fusion relations for 7(z) imply the following recursive relations for Ti(z): 

Tl(z)n(zq) = T£ + l(z) + T£„l(z(J 2 ) 

which can be solved in terms of determinants [5] : 

/ n(z) 1 
1 n(zq) '•• 



rt(z) = det 



1 



V o 



l n(V _1 ) / 



The remarkable fact is that elements {rg(z)} also satisfy another set of relations 
which also follow from the fusion relations: 



Ti(zq*)ri(zq 2 ) = n + i(z)n-i{z) + 1. 



(28) 



18 



3.7 Local integrable quantum spin Hamiltonians 

Transfer-matrices 

t m (u) = MKr^/ai) ■ • • icrwai)4°) 

form a commuting family of operators in y( mi ) <g) • • • <g) y( m «). They quantize 
the generating functions for classical spin chains and can be used to construct 
local quantum spin chains. Below we outline two common constructions of local 
Hamiltonians from such transfer-matrices. 

3.7.1 Homogeneous SU(2) spin chains 

The homogeneous Heisenberg model of spin S corresponds to the choice mi = 
• ■ • = mjv = I and = ■•• = &i = 1- We will denote corresponding transfer 
matrices as (u) 

In the case m = 1, the transfer matrix t^\z) is 

t { l\z) =tro(RoN(z)...Roi(z)) (29) 

where R(z) is the matrix (fT9"|) . 

The linear operator (|29|) is the transfer matrix of the 6-vertex model [23] [3] . 
It is also a generating function for local spin Hamiltonians: 

, n 
H^—hg^iu))^ = £(a>? 1+1 +« +1 + A<< +1 ), (30) 

ra=l 

, N 
H k = (£)* logfW^Uo = # ( *V»> • ■ • ,Vn+k). (31) 

n— 1 

Here a x ,a y ,a z are Pauli matrices and a n is the collection of Pauli matrices 
acting nontrivially in the n-th factor of the tensor product. 

-■-(;;)■ H -.o)< ---(j-i)- <*> 

If l > 1 similar analysis can be done for the transfer-matrix 

tf\z) = tr a (R^\z)...R«$(z)) 

Since R^ l \l) = P, we have: 

tf\l) = tr a (P al . . . P aN ) = tr a (P 12 P 13 . . . P 1N P al ) = P 12 P 13 ...Pin 

Here we used the identities P a \A a P a \ — Ai and tr a (P a i) = 1^. 
The operator T — t\ (1) is the translation matrix: 

T(xi (g> X2 ■ ■ ■ <8> x n ) = xn <8> xi <E) ■ ■ ■ £jv-i 
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Differentiating tf\z) at z — 1 we have: 



N N 



Similarly, taking higher logarithmic derivatives of tj (z) at z = 1 we will 



have higher local Hamiltonians acting in (C i+1 ) 5 



Here the matrix acts in (C' +1 )® fc+1 . The subindices show how this matrix 
acts in (C m )® jv . 

One can show that these local quantum spin chain Hamiltonians in the limit 
q — > 1 and I — > co become classical Hamiltonians described in section 12.41 
assuming that q l is fixed. 

3.7.2 Inhomogeneous SU(2) spin chains 

The construction using degeneration points. 

The construction of inhomogeneous local operators is easy to illustrate on 
the spin chain where the inhomogeneities alternate. 

t m (z) = tr a (R^ M (za-i)R% M (za) . . . R^Uza^R^za)) 
Now we have two sublattices and two translation operators 

Teven = ^24^26 ■ • • P2,2N, T a dd = fl3-Pl5 ■ • • P\,2N-l 

It is easy to find the following special values of the transfer-matrix: 
t h (a) = T™™Rt h) (a- 2 ) . fl^-itO 

UAa- 1 ) = Rt h \a 2 ) ■ ..Rt% N (* 2 )T° dd 
These operators commute and 

*ii(o)*fa(o _1 ) =T even T odd 
t^a^ia- 1 )- 1 =T even (T° dd )- 1 

Tf{h,h) ( n -2\T>(l2,h){-2\ r>(h,h) , -2\ r,(h ,h) („-2\ u(h,h) f„-2\ 
""Z^N—IV* ^-"-4,1 \ u I ■■■ ■ ft 2iV,2./V-3 \ u ) JX 2N-1,2N\"' ) • ■ ■ ^1,2 \ u I 

(33) 



20 



Taking logarithmic derivatives of (z) at z — a and of ti 2 (z) at z = a 1 we 
again will have local operators, for example: 



These Hamiltonians in the semiclassical limit reproduce inhomogeneous clas- 
sical spin chains described earlier. 

4 The spectrum of transfer-matrices 
4.1 Diagonalizability of transfer- matrices 

Assume that q is real. Let t(u)* be the Hermitian conjugation of t(u) with 
respect to the standard Hermitian scalar product on (C 2 )®^. It is easy to 
prove, using the identities for R(u) that 



where z is the complex conjugate to z. 

Because t(z) is the commutative family of operators, the operators t(z) is 
normal when ~ai — 1 . Therefore for these values of a, it is diagonalizable. Since 
t(z) is linear (up to a scalar factor) in af, this imply that t(z) is diagonalizable 
for all generic complex values of di, and for the same reasons for all generic 
complex values of q. 

4.2 Bethe ansatz for sl 2 

In this section we will recall the algebraic Bethe ansatz for the inhomogeneous 
finite dimensional spin chain. 

The quantum monodromy matrix for such spin chain is: 



Z-T-l0gt h {z)\ z = a = ^2R 2 ^ + l,2n( a2 ) R '2n+l,2n( a2 ) 




t(z\a u ...,a N )* = {-l^tiz^q- 1 ^- 1 . . .aW' 1 ) 



T(z)=L ( i mi \z/a 1 )...C^ N \z/a N )D 



(35) 



where 




It is convenient to write it as 
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In the basis e\ ® e\, e\ ® e2, e-i ® e.\, ® &2 of the tensor product C 2 ® C 2 wc 
have: 



Ti(zw)T 2 (w) 



/ A(zw)A(w) A(zw)B(w) B(zw)A(w) B(zw)B{w) \ 

A(zw)C(w) A(zw)D(w) B(zw)C(w) B(zw)D(w) 

C{zw)A{w) C(zw)B(w) D(zw)A(w) D(zw)B(w) 

\ C{zw)C{w) C(zw)D(w) D(zw)C(w) D(zw)D(w) J 



Writing the i?-matrix in the tensor product basis as in (fl 



/ 1 \ 

f{z) g(z)z- 1 

g(z)z f(z) 

\ 1 / 



the commutation relations (1161) produce the following relations between A and 
B and D and B: 



A(z)B(v) 
D(z)B(v) 



l -,B{v)A {z) - B {Z) A( V ) 



f(vz- r ) 

1 



f(zv 



B{v)D{z)- 9 -^^B{:)D [ r ] 



where f(z) 



-i ,g(z) 

zq—z 1 q 1 ' ^ v / 



— 1~1 



The L-operators act on the vector 



i mi) ®v ( m2) 



fizv- 1 ) 



(ton) 
^0 



in a special way: 



Li(z)£l 



{zq- 



z x q 




(zq~ 



z 1 q 



From here it is clear that fi is an eigenvector for operators A and D and 
that C annihilates it: 



T(z)n 



a(z)fl * 
S(z)Q 



where 



N 



a(z) = zY\_(za { x q~ 



z 1 a i q 



i=l 
N 



6(z) = Z- 1 H(za; 1 c 



m i - 1 

z~ L aiq^~ ) 



The details of proof of the following construction of eigenvectors can be 
found in IT31. 
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Theorem 1. The following identity holds 

(A(z) + D(z))B( Vl ) . . . B(v n )Sl = A(2|{«,})B(«i) . . . B(v n )Sl 

where 

«iZ g — V: zq -f-j v i zq - ViZ q 



a / I r / \ TT iZ 1~ v i z 1 , r/ ^TT t, ^ zq- Vi z q 

Hmvi}) = a{z) [[ — +S(z)\\^—— L — (36) 

if the numbers Vi satisfy the Bethe equations: 



TT ~ ^ = _ z2 tt ^ ;g - "j 'wr (37) 

Note that the formula for the eigenvalues in terms of solutions to Bethe 
equations is a rational function. Bethe equations can be regarded as conditions 

res z=Vj A(z\{vi}) = 

This agrees with the fact that t(z) is a commuting family of operators which 
has no poles at finite z. 



4.3 The completeness of Bethe vectors 

The next step is to establish whether the construction outlined above all eigen- 
vectors. We will focus here on the spin chain of spin 1/2. 

Assume that q, e 2H , and inhomogeneities are generic. Let us demonstrate 
that the vectors 

B(v 1 )...B(v n )Q (38) 

where Vi are solutions to Bethe equations give all 2 N eigenvectors of the transfer- 
matrix. 



4.3.1 

Consider the limit of (|38l) when ajv — > oo. 

Assume that v is fixed and —¥ oo. From the definition of B(v) we have: 

B(v) = a N q~^v~ l (A{v) ® / - B(v) ® K){1 + o(l)) 

where A(v), B{v) are elements of the quantum monodromy matrix (|35p with 
only N — 1 first factors. 

On the other hand if — ¥ oo and v — > oo such that v — wa^ and w is 
finite the asymptotic is different: 

N 

B(v) = w~ 1 q~i Y[(wa N a~ l qi)k ® • • • ® ® /(l + o(l)) 

n=l 
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From here we obtain the asymptotic of Bethe vectors when all Vi are fixed 
and cin — > oo 



B( Vl ) . . . B(v n )n N -> q^ n (-a N qi) n J] 

n 

(B( Vl ) . . . B(v n )n N ^ ® - q-^- 1 ' 1 

=1 



B(vi) . . . A(vi) . . . B(v n )tl N -i ® /^ mw) )(l + o(l)). (39) 

Similarly, when v±, . . . ,v n —i are fixed and a at — ¥ oo such that v n = wa^ we 
have 

n-l 

B(vi) . . . B(v n )Q N -> ^^^-"(-W)"-! J] ,,-lgT 1 

i=l 

B(«0 . . . B(v n -i)Sl N -i ® /^ mjv) (l + o(l)) (40) 

4.3.2 

Solutions to (|37p have the following possible asymptotic when ajy — > oo : 

1. For all j = 1,...,N, lim aN ^ 00 Vj — v'a where {v'j} is a solution to 
the Bethe system for the spin chain of length N — 1 with inhomogeneities 
ai, . . . , cln-i and Z. 

2. For one of Vj's, say for v n we have v n = cinw + 0(1) and for others 
lim ajv ^oo = v'j where {v'A is a solution to the Bethe system for the spin 
chain of length N — 1 with inhomogeneities oi, . . . , a^-i and Zq^ 1 . From the 
Bethe equation for v n we have 



w 2 



1 - Z 2 q- N+2n 
q 2 - Zq- N + 2n 



3. More then one of Vi is proportional to a at. 

Using induction and the asymptotic of Bethe vectors (|39| and (|40]) it is 
easy to show that only first two options describe the spectrum of the spin 1/2 
transfer-matrix. Similar arguments were used in |18j to prove the completeness 
of Bethe vectors in an SL n spin chain. 

( N \ 

This implies immediately that there are I J Bethe vectors for each < 
n < N. And that the total number of Bethe vectors is 

N 



2" = E 



n=0 



N 
n 



Other solutions to Bethe equations describe eigenvectors in infinite dimensional 
representations of quantized affine algebra with the same weights. They do not 
correspond to any eigenvectors of the inhomogeneous spin 1/2. 



24 



4.3.3 



For special values of a a the solutions to the Bethe equations may degenerate (a 
level crossing may occur in the spectrum of t(z)). In this case the Bethe ansatz 
should involve derivatives of vectors (|3"51) . 

5 The t her mo dynamical limit 

The procedure of " filling Dirac seas" is a way to construct physical vacua and 
the eigenvalues of quantum integrals of motion in intcgrable spin chains solvable 
by Bethe ansatz. 

To be specific, consider the homogenous spin chain of spins 1/2. Let Hi, H2, ■ ■ ■ 
be quasilocal Hamiltonians described in the section 13.7.11 with q = exp(r?) for 
some real r\. 

Take the linear combination 

H(X)^Y, H ^ ( 41 ) 

k 

This operator is bounded. Let fijv(A) be its normalized ground state. As 
N — > 00, matrix elements (Q,n, ail at) converges to the state lu\ on the inductive 
limit of the algebra of observables. The action of local operators on lu\ generate 
the Hilbert space "K. 

Since the eigenvalues of coefficients of t(u) can be computed in terms of 
solutions to Bethe equations, the spectrum of these operators in the large N is 
determined by the large N asymptotic of solutions to Bethe equations. 

The main assumption in the analysis of the Bethe equations in the limit 
N — > 00 is that the numbers {vi°^}0are distributed along the real line with some 
density p(u). The intervals where p(u) ^ are called Dirac seas. For Hermitian 
Hamiltonian (|4"T|) there is strong evidence that Dirac seas is a finite colelction of 
intervals (B^~ , -Bf), . . . , (-B+, B~). Here numbers B^ are boundaries of Dirac 
seas. We assume they are increasing from left to right. The boundaries of Dirac 
seas are uniquely determined by {A/}. 

A solution to the Bethe equations is said to contain an m-string, when as 
N — >• 00, there is a subset of {vi} of the form 

v (m) +i%n 1 v (m) +i^{m- 2),...y m ) - i|(m - 2), v (m) - i~m. 

with some real v^ n K 

The excitations over the ground states can be of the following types: 

• A hole in the Dirac sea (B% , B^ ) correspond to the solution to Bethe 
equations which has one less number Vi, and as N —¥ 00 the remaining 
Vi "fill" the same Dirac seas with the the densities deformed by the fact 
that one of the numbers {v^} is missing and other are "deformed" by 

Solutions to the Bethe equations corresponding to the minimum eigenvalue of H(X) 
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the missing one. The number which is "missing" is a state with a hole is 
the "rapidity" v € (B%,BZ) of the hole. 

• Particles correspond to "adding" one real number to the collection 

• m-strings, m > 1 (corresponding to adding one m-string solution to the 
collection {va ^})- 

There are convincing arguments, that the Fock space of the system with 
the Hamiltonian (|41 [) has the follwoing structure. It has a vacuum state 0\ 
corresponding to the solution of the Bethe equations with the minimal eigenvaue 
of H(X). Excited states are eigenvectors of the Hamiltonian (and of all other 
integrals) which corresond to solutions of Bethe euqations with finitely many 
holes, particles, and m-strings. It has the following structure 

^({5+ sail) =0 

A' h >0,A'p>0,Ar>0 ni + ---+n k =N h ,N p 

■ (g) Lf mm (/+ xni x . . . /+ xnk x I N * x s lN ) 

Here we used the notation 7 ; + = (£? ; + , B7), 17 = {Bf_,, B7), where Bq = —tt, 
and B7, 1 = n; is the number of holes in the Dirac sea {B+,B~), N p is the 
number of particles; I is the complement to the Dirac seas, N is the number 
of strings with m > 2. The symbol "symm" means certain symmetrization 
procedure which we will not discuss here (see, for example, |17| for a discussion 
of the aniferromagnetic ground state) . 

Varying {Afc}, or equivalently, positions {B k } of Dirac seas we obtain a 
"large" part of the space of states of the spin chain in the limit N —> oo: 
(C 2 )®^ in the limit N —> oo into the direct integral of separable Hilbert spaces: 

(C 2 )® Ar ^0r K{{Bt,BT}U) (42) 

6 The 6-vertex model 

6.1 The 6-vertex configurations and boundary conditions 

The 6-vertex model is a model is statistical mechanics where states are con- 
figurations of arrows on a square planar grid, see an example on Fig. [5] The 
weights are assigned to vertices of the grid. They depend on the arrows on 
edges surrounding the vertex non-zero weights correspond to the configurations 
on Fig. [TJ to configurations where the number of incoming arrows in equal to 
the number of outgoing arrows. This is also known as the ice rule [24] . 

Each configuration of arrows on the lattice can be equivalently described 
as the configuration of "thin" and "thick" edges (or "empty" and "occupied" 
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a, a 2 b! b 2 c } 




Figure 1: The 6 types of vertices and the corresponding thin and thick edges 
configurations. 

edges) as it is shown on figure [TJ There should be an even number of thick edges 
at each vertex as a consequence of the ice rule. 

The thick edges form paths going from NorthWest (NW) to SouthEast (SE). 
We assume that when there are 4 think edges meeting at a vertex, the corre- 
sponding paths meet at this point and then are going apart. So, equivalently, 
configurations of the 6- vertex model can be regarded as configurations of paths 
going from NW to SE satisfying the rules from Fig. Q0. 

6.2 Boundary conditions 

It is natural to consider the 6- vertex model on surface grids. 

If the surface is a domain on a plane we will say that an edge is outer if it 
intersects the boundary of a domain. We assume that the boundary is chosen 
such that it intersects each edge at most once. Outer edges are attached to a 
4-valent vertex by one side and to the boundary by the other side. 

Fixed boundary conditions means that fixing the 6- vertex configurations on 
outer edges. An example of fixed boundary conditions known as domain wall 
(DW) boundary conditions on a square domain is shown on fig. [5] 

We will be interested in three types of boundary conditions: 

• A domain (connected simply connected on a plane) with fixed boundary 
conditions, see Fig. [2j We will also call this Dirichlet boundary conditions. 

• A cylinder with fixed boundary conditions, see Fig. ??. This case can 

One can consider such configurations on any 4-valent graph. But only for special graphs 
and special Boltzmann weights one can compute the partition function pe site. 
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Figure 2: A domain, connected, simply-connected. 



Figure 3: Cylindric boundary conditions. 

be regarded as a domain with states on outer edges of two sides being 
identified and with fixed boundary conditions on other sides. 

• Identification of states on outer edges of opposite sides of a rectangle gives 
the states for the 6-vertex model on a torus, see Fif. [4] It is also known 
as the 6-vertex model with periodical boundary conditions. 

6.3 The partition function and local correlation functions 

To each configuration oi, 02, 61, 62, ci, and C2 on Fig. [T] we assign Boltz- 
mann weights which we denote by the same letters. The physical meaning of a 
Boltzmann weight is exp(— ^), where E is the energy of a state and T is the 
temperature (in appropriate units), so all numbers a\, 0,2, b\, &2, Ci, and ci 
should be positive. 

The weight of the configuration is the product of weights corresponding to 
vertices inside the domain of weights assigned to each vertex by the 6-vertex 
rules. 

The 6-vertex model is called homogeneous if the weight assigned to a vertex 
depends only on the configuration of arrows on adjacent edges and not on the 
vertex itself. 

When the weights also depend on the vertex the model is called inhomoge- 
neous. 
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Figure 4: Toric boundary conditions. 



The partition function is the sum of weights of all states of the model 
Z = JJ w (vertex), 



states vertices 



where iu(vertex) is the weights of a vertex assigned according to Fig. [TJ 

Weights of states define the probabilistic measure on the set of states of the 
6-vertex model. The probability of a state is given by the ratio of the weight of 
the state to the partition function of the model 

P{state) = nvertiees ^(vertex) ^ ^ 

The characteristic function of an edge e is the function defined on the set of 
6-vertex states 



<r e (state) 



1, e is occupied by a path, 
0, otherwise. 



A local correlation function is the expectation value of the product of such 
characteristic functions with respect to the measure (|43[) : 

n 

(a ei a e2 ..a en ) = P(state) a ei (state). 

states i—1 

It is convenient to write the Boltzmann weights in exponential form 

H+V -H-V 

<2i = ae , a2 = ae , 
b x = be+ H ~ v : b 2 = be- H+v , 

-E E 

c\ = ce , C2 = ce , 

From now on we will assume that E = 0. If the weight are homogenous (do 
not depend on a vertex), local correlation functions for a domain, cylinder or 
torus do not depend on E. Also, the parameters H and V have a clear physical 
meaning, they can be regarded as horizontal and vertical electric fields. Indeed, 
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Figure 5: A possible configuration of paths on a 5 x 5 square grid for the DW 
boundary conditions. 

if E = the weight of the state can be written as 
w(S) = Yl w(v\S) 

v&vertices 

exp( £ a e (S)(H(e + ) + ff(e_))/2 + £ cr e (S)(V(e+) + V(e_))/2) (44) 

e£E h e£E v 

Here S is a state of the model, Eh is a set of horizontal edges and E v is the 
set of vertical edges, w(v\S) is the weight of the vertex v in the state S where 
a\ = a-2 = a and b\ = bj = b. The symbol o~ e (S) is the characteristic function 
of e: <J e {S) — 1 if the arrow pointing up or to the left, and <r e {S) = —1 if the 
arrow is pointing down or to the right. 

For a given domain let us chose its boundary edge and enumerate all other 
boundary edges counter-clock wise. The partition function for a domain for 
various fixed boundary conditions can be considered as a vector in (C 2 )® N 
where the factors in the tensor product counted from left to right correspond to 
enumerated boundary vertices. 

Similarly, the partition function for a (vertical) cylinder of size N x M (hor- 
izontal, vertical) can be regarded as a linear operators acting in (C 2 ) 87V where 
the factors in the tensor product correspond to states on boundary edges. 

The partition function for the torus of size N x M is a number and it can be 
regarded as a the trace of the partition for the vertical cylinder of size N x M 
over the states on its horizontal sides, i.e. over (C 2 )® N . It can also be regarded 
as trace of the partition for the horizontal cylinder of size M x N over the states 
on its vertical sides, i.e. over (C 2 )® M . The result is an identity which we will 
discuss later. 
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6.4 Transfer-matrices 



Let is write the matrix of Boltzmann weights for a vertex as the 4x4 matrix 
acting in the tensor product of spaces of states on adjacent edges. 

Let e\ be the vector corresponding to the arrow pointing up on a vertical 
edge, and left on a horizontal edge. Let e 2 be the vector corresponding to the 
arrow pointing down on a vertical edge and right on a horizontal edge. 

The matrix of Boltzmann weights with zero electric fields acts as 



Rei <x> e\ 
Rei ® e 2 
Re 2 <8> ei 
Re 2 <8> e 2 



ae\ (g> ei 

bei <£> e 2 + ce 2 ® e\ 
be 2 <8> ei + cei ® e 2 
ae 2 ® e 2 



In the tensor product basis e\ ® e\, e\ ® e 2 , e 2 <8> ei, e 2 i 



/ 


o 











\ 









c 












c 


6 







V 











a 


J 



(45) 
(46) 
(47) 
(48) 

e 2 it is the 4x4 matrix 
(49) 



The 6- vertex rules imply that the operator R commutes with the operator rep- 
resenting the total number of thick vertical edges, i.e. 

[D H <E>D H ,R] 







where 



D 



H 



e H/2 







-H/2 



The row-to-row transfer-matrices with open boundary conditions also known 
as the (quantum) monodromy matrix is defined as 

T - D H{aS) B D H(a.l) jjH(a,l) r> 

It acts in the tensor product <£ 2 a ® Cf ® • • • <£> C^r of spaces corresponding to 
horizontal edges and vertical edges . Each matrix R a i is of the form (??), it 
acts trivially (as the identity matrix) in all factors of the tensor product except 
a and i. In the inhomogeneous case matrix elements of R a i depend on i. 

A matrix element of T is the partition function of the 6- vertex model on a 
single row with fixed boundary conditions. 

Define operators 

D (a) = D V(a,l) . . . D V(a,N) 



D" = l 



D 



H 



1 



The row-to-row transfer-matrix corresponding to the cylinder with a single 
horizontal row a with with electric field H (a, i) applied to the i-th horizontal 
edge of the a-th horizontal line is the following trace 

t a = tr a (T a ) = tr(D«^R al . . . R aN ) 
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It is an operator acting in (<C 2 )® N . Its matrix element is the partition function 
of the 6-vertex model on the cylinder with a single row with fixed boundary 
conditions on vertical edges. 

The partition function of an inhomogeneous 6-vertex on a cylinder of height 
M with fixed boundary condition on outgoing vertical edges is a matrix element 
of the linear operator 

where 

D (a) = jjV(a.l) ^ . . . D V{a,N) 

The partition function for the torus with N columns and M rows is the trace 
of the partition function for the cylinder 

Z nm = tr iC 2 )e>K (D^ M h M ■ ..D^h) 
Using the 6-vertex rules this trace can be transformed to 

Znm = tr {c *)®»(t { M M) • ••*i ifl) (£>i /l ® ••• ® D%*)) 

where H a = YliLi H ai and = £0=1 Ki and 

=tr a {R aX ...R aN D*) 

The partition function for a generic domain does not have a natural expres- 
sion in terms of a transfer-matrix. However, it is possible in few exceptional 
cases, such as domain wall boundary conditions on a square domain, sec for 
example [16] . and references therein. 



6.5 The commutativity of transfer-matrices and positivity 
of weights 

6.5.1 Commutativity of transfer-matrices 

Baxter discovered that matrices of the form (??) acting in the tensor product 
of two two-dimensional spaces satisfy the equation 

R12R13R23 = #23-^13-^12 (50) 

if 

a 2 + b 2 - c 2 a' 2 + b' 2 - c' 2 a" 2 + b" 2 - c" 2 



2ab 2a' V 2a"b" 

This parameter is denoted by A: 

, a 2 + b 2 - c 2 



2ab 

If each factor in monodromy matrices T' a — R' al . . . R' aN and T' h ' = . . . R'^ N 
have the same value of A, the equation (|50l) implies that they these monodromy 
matrices satisfy the relation 

n ab l a l b - l b l a H ab 
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in V a ® Vb ® Vi ® ■ ■ ■ (g) Vat. 

If R is invertible, this relation implies that row-to-row transfer-matrices with 
periodic boundary conditions commute: 

t' =tr a {T' a D*), t" = tr b (Tl'D*), [t',t"] = 

It is easy to recognize that t is exactly the generating function for commut- 
ing family of local Hamiltonians for spin chains constructed earlier. Thus, the 
problem of computing the partition function for periodic and cylindrical bound- 
ary conditions for the 6- vertex model is closely related to finding the spectrum 
of Hamiltonians for integrable spin chains. 

6.5.2 The parametrization 

The set of positive triples of real numbers a : b : c (up to a common multiplier) 
with fixed values of A has the following parametrization [3]. 

1. When A > 1 there are two cases. 

If a > b + c, the Boltzmann weights a, 6, and c can be parameterized as 

a = r sinh(A + 77), 6 = rsinh(A), c = rsmh(rj) 
with A, rj > 0. 

If b > a + c, the Boltzmann weights can be parameterized as 

a = r sinh(A — 7/), 6 = rsinh(A), c = rsinh(?7) 

with < r\ < A. For both of these parameterizations of weights A = 
cosh(?7). 

2. When -1< A < 

a = rsin(A — 7), 6 = rsin(A), c = rsin(7), 
where < 7 < 7r/2, 7 < A < tt/2, and A = — COS7. 

3. When < A < 1 

a = rsin(7— A), 6 = rsin(A), c = rsin(7), 
where < 7 < 7r/2, < A < 7, and A = cos 7. 

4. When A < — 1 the parametrization is 

a — r sinh(ry — A), fo = rsinh(A), c = rsinh(?7), 
where < A < r/ and A = — cosh 77. 
We will write a = a(w), b = b(u), c — c{u) assuming this parameterizations. 
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Figure 6: The Yang-Baxter equation. 




Figure 7: The 'unitarirty'. 

6.5.3 Topological nature of the partition function of the 6-vertex 
model 

Fix a domain and a collection of simple non-selfintcrsccting oriented curves with 
simple (transversive) intersections. The result is a 4-valent graph embedded into 
the domain. Fix 6-vertex states (arrows) at the boundary edges of this graph. 
Such graph connects boundary points, and defines a perfect matching between 
boundary points. 

The 6-vertex rules define the partition function on such graph. The Yang- 
Baxter relation implies the invariance with respect to the two moves shown of 
Fig. El Fig. [3 Because of this, the partition function of the 6-vertex model de- 
pends only on the connection pattern between boundary points, i.e. it depends 
only on the perfect matching on boundary points induced by the graph. 

The 'unitarity' relation involves the inverse to R, and therefore does not 
preserve the positivity of weights. But if it is used "even" number of times it 
gives the equivalence of partition functions with positive weights. For example 
it can be used to 'permute' the rows in case of cylindrical boundary conditions. 

6.5.4 Inhomogeneous models with commuting transfer-matrices 

From now on we will focus on the 6-vertex model in constant electric fields. 
If Boltzmann weights of the 6-vertex model have special inhomogeneity ay — 
a(v,i — Wj), bij = b(ui — Wj), otj = c(m — Wj) the partition function of the 6-vertex 
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model on the torus is 

4 T l(M> Ml#, *0 = (*(«i) ■ • • i(«M)(fl W ® ■ ■ ■ ® (51) 

where Wjv - (C 2 )® w and i(u) - tr a (T a (u)) 

T a {u) = R al {u - Wl ) . . . R aN {u - w N )D NH (52) 
Because the i?-matrix satisfies the Yang-Baxter equation we have 
R ab (u - v)T a (u)T b (v) = T b (v)T a (u)R ab (u - v) 
As a corollary, traces of these matrices commute: 

[t(u),t(«)] =0 

and as we have seen the previous sections their spectrum can be described 
explicitly by the Bcthc ansatz. 

Notice that positivity of weights restricts possible value of inhomogeneities. 
For example when A > 1 and a > b + c we should have —u<Wi<u. 

6.5.5 

For a diagonal matrix d such that [d ® 1 + 1 ® d, R(u)] = define 
R d (u) = (e ud ® l)R(u)(e- ud <g> 1) 

It is clear that 

tr(R d al (u- Wl ) . . . R d aN {u-v N )D NH ) = Utr{R al {u- Wl ) . . . R aN (u-v N )D NH )U- 1 

where U = exp(- Y^,iLi w i d i)- 

In particular the partition function for a torus with weights given by R(u) 
and with weights given by R d (u) are the same. 

6.5.6 

Let u\,...,Un be parameters of inhomogeneities along horizontal directions, 
and wi, . . . ,Wm be inhomogeneities along vertical directions. The partition 
function (??) have the following properties: 

• Z N M ({u}, {w}\H, V) is a symmetric function of u and w. 

• It has a form njl e ~ Mu * Utii e~ Nwi P({e 2u }, {e 2w }) where P is a poly- 
nomial of degree M in each e Ui and of degree TV is each e Wj . 

• It is a function of m — Wj . 

• It satisfies the identity 

4 T l({«l, {w}\H, V) = Z$ N ({w}, {h - u}\V, H) (53) 

where h — r\ or 7 depending on the regime. This identity correspond to the 
"rotation" of the torus by 90 degrees and is known as "crossing-symmetry" 
or " modular identity" . 
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7 The 6- vertex model on a torus in the thermo- 
dynamic limit 

7.1 The thermodynamic limit of the 6- vertex model for 
the periodic boundary conditions 

By the thermodynamical limit here we will mean here the large volume limit, 
when N, M — > oo. 

The free energy per site in this limit is 

log(Z N , M ) 

N,M^-oo NM 

where Zn^m is the partition function with the periodic boundary conditions 
on the rectangular grid Ln,m- It is a function of the Boltzmann weights and 
magnetic fields. 

For generic H and V the 6- vertex model in the thermodynamic limit has the 
unique translational invariant Gibbs measure with the slope (h,v): 

1 df 1 1 df 1 

h = — TT- + -, v= — + (54) 
2dH 2' 2dV 2 V ; 

This Gibbs measure defines local correlation functions in the thermodynam- 
ical limit The parameter 

a 2 + b 2 - c 2 



A 



2ab 

defines many characteristics of the 6- vertex model in the thermodynamic limit. 

7.2 The large N limit of the eigenvalues of the transfer- 
matrix 

The row-to-row transfer-matrix for the homogeneous 6-vertex model on a lattice 
with periodic boundary condition with rows of length N is 

N 

t(u) = tr a {R al {u) . . . R aN (u)e H ^)exp{V ^ a z a ) 

a=l 

According to (|3"6"|) and (|3"T1) . the eigenvalues of this linear operator are: 

A { „- } = a N e NH IT VlZ ~ lq Vil "f l + b*e- NH TT V ^ Z \ - (55) 
f = i Viz - Vi z f = \ Vi z - Viz 

where < n < N and z,q are parameterizing a, b, c as a = r(zq — zq^ 1 ,b — 
r(z — z~ x ),c — r(g — The numbers Vi are solutions to the Bethe equations: 



'•«-"-'<r\ = _ e - M fr!!iWy (56) 

~=i ViVj q 1 - v t vjq 



Vi - v, 1 
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As for the corresponding spin chains it is expected that the numbers Vi 
corresponding to the largest eigenvalue concentrate, when N — > oo, on a contour 
in a complex plane with the finite density. The Bethe equations provide a linear 
integral equation for this density. This conjecture is supported by the numerical 
evidence and it is proven in some special cases, for example when A = 0. 

The partition function for the homogeneous 6-vertex model on an TV x M 
lattice with periodic boundary conditions is 

a 

where a parameterize eigenvalues of t(u). 

Let ojtv be the eigenvector of t(u) corresponding to the maximal eigenvalue. 
The sequence of vectors {VL N } as N — > oo defines the Hilbert space of pure states 
for the infinite system. Let A be the largest eigenvalue of t(u). According to 
the main conjecture about that the largest eigenvalue correspond to numbers 
Vi filling a contour in a complex plane as JV, the largest eigenvalue has the 
following asymptotic: 

AW=exp(-iV7(#,F) + 0(l)) (58) 

The function f(H, V) as we will see below it is the free energy of the system. It 
is computed in the next section. 

The transfer-matrix in this limit has the asymptotic 

t(u) = exp(-Nf(H,V))t(u) 

where the operator t(u) acts in the space and its eigenvalues are determined 
by positions of "particles" and "holes", similarly to the structure of excitations 
in the large N limit in spin chains. 

7.3 Modularity 
7.3.1 

It is easy to compute now the asymptotic of the partition function in the ther- 
modynamical limit. As M — > oo the leading term in the formula (|57|) is given 
by the largest eigenvalue: 

Z N<M =A^ M (l + 0(e- M )) 

for some positive a. 

Taking the limit N — > oo and taking into account the asymptotic of the 
largest eigenvalue (|58|) we identify the function f(H, V) in (|58j) with the free 
energy: 

z NM = e ™/ew<i+°ci)) 
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7.3.2 



Notice now that we could have changed the role of N and M by first taking the 
limit N — > oo and then M — > oo. 

In this case we would have to compute first the asymptotic of 

Z n ,m = tr(t(u) M ) 

as N — > oo and then take the limit M —> oo. 

The large N limit of the trace can be computed by using the finite temper- 
ature technique developed by Yang and Yang in It was done by de Vega 
and Dcstry [10 . The leading term of the asymptotic can be expressed in terms 
of the solution to a non-linear integral equation. 

This gives an alternative description for the largest eigenvalue. Similar de- 
scription exists for all eigenvalues. In other integrable quantum field theories it 
was done by Al. Zamolodchikov [40 . 

8 The 6-vertex model at the free fermionic point 

When A = the partition function of the 6-vertex model can be expressed in 
terms of the dimer model on a decorated square lattice. Because the dimer 
model can be regarded as a theory of free fermions, the 6-vertex model is said 
to be free fermionic when A = 0. 

At this point the raw-to-raw transfer-matrix on TV-sites for a torus can be 
written in terms of the Clifforrd algebra of C^. The Jordan- Wiegner transform 
maps local spin operators to the elements of the Clifford algebra. 

In Bcthe equations (|37[) the variables Vi disappear in the r.h.d which be- 
comes simply (—l) n ~ 1 e~ 2NH . After change of variables these equations can be 
interpreted as the periodic boundary condition for a fermionic wave function. 

8.1 Homogeneous case 
8.1.1 

At the free fermionic point A = 0, i.e. 7 = § and the weights of the 6-vertex 
model are parameterized as: 

7T 

a = sin(— — u ) = cos(u), b — sin(u), c = 1 

Without loosing generality we may assume that < u < ? . 

The eigenvalues of the row-to-row transfer-matrix are given by the Bethc 
ansatz formulae: 

n 

A(u) = ((cos{u) N (-l) n e NH + {sm{u)) N e- NH )e NV J^(cot(u - v 4 )e~ 2V ) (59) 

»=i 
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Here < n < N, and Vi are distinct solutions to the Bethe equations: 

cot(v) N = (-lf-i e - 2Jvff 



or: 



cot(w) = uie 



-2H 



where uj n = 

Using the identity: 



cot u cot v + 1 1 — cot 2 u 

cot(u — v) = = cot u — 



cot v — cot u cue 2H — cot u 

we can write the eigenvalues as: 

A(u) = ((cos( U )^(-l)" e Nff +(sin( U )) Ar e - Jvff ) e ( Ar ~ 2 ») l/ )fT(cot U 1 ^ - ) 

uJie~ ZM — cot u 

i — 1 



(60) 



8.1.2 



To find the maximal eigenvalue in the limit N —¥ oo we should analyze factors 
in the formula (|59p . 

1. If 

max cot(w — v) < e 2V 

cj| — 1 

where cot(u) = we _2fl , all factor are less then one and the maximal eigenvalue 

A ord (u) = (cos(u) N e NH + (sm(u)) N e- NH )e NV 

is achieved when n = 0. 

As N — > oo the first term dominates when cot(u) < e 2H . In this case the 
Gibbs state describing the 6- vertex model in the thermodynamical limit is the 
ordered state A\ from Fig. ITT1 When cot(w) > e 2H , the second term dominates. 
In this case the Gibbs state the ordered state B-i shown on Fig. [TT] 

2. If mini^^i cot(ii — v) > e 2V all factors cot(u — Vj)e 2V are greater then 
one by absolute value. In this case the maximal eigenvalue is achieved when 
n = N. The corresponding ground state is ordered and in A2 from Fig. [Til 
when cotw > e~ 2H and is B\ when cot(u) < e~ 2H 

3. If 

max cot(w — v) < e 2V 

cj| — 1 

then there exists u>q — e lK such that 

I cot(u — 6)1 =e 2y 
where cot(6) = e ± lK - 2H . It is easy to find K: 

(Ue 2H + U- 1 e- 2H )e iV - (Ue- 2H + U^e 211 ) 



cos(K) = 



2(e 4V + 1) 
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In this case 



| cot(u — v)\ > e 



2V 



when cot(u) = e la e with —K < a < K and 



2V 



| cot(u — v)\ < e 

when — 7r < a < —K or K < a < 7r. 

The maximal eigenvalue is this case corresponds to maximal n such that 



7r(n— 1) 



— < K and w 3 = exp( 



■ 7r(n+l-2jh „■ _ 



iV 



) , j = 1 , . . . , n and is given by §5 



As N — > oo the asymptotic of the largest eigenvalue is given by the integral: 



\og(A disord (uj) = A(log(cos(u)e ff ) + 



1 

2ttI 



K 



log 



/ cot(u)we 2H + 1\ duj 



-K 



\ ue — cot(it) / u 



-)+o(i) 



The 6-vertex in this regime is in the disordered phase. 

Disordered and ordered phases are separated by the curve 



max | cot(it — v)\ = e , min | cot(it — v)\ = e 
M=i | W |=1| 



-IV 



or, more explicitly 



Ue 



-2H 



U ± e~ 2H 



,-2V 



This curve is shown on Fig. [8] 



8.2 Horizontally inhomogeneous case 

The partition function for the 6-vertex model at the free fermionic point with 
inhomogeneous rows is 



Zn,m( u i a ) = ^ A(u + a\uj) M A(u - a\u) 

n— u;i,...,o; n 



M 



where w,; 7^ cjj are solutions to 



(-1) 



n-l 



and A(u|w) is given by (|6"U|) . 

The boundary between ordered and disordered phases is given by equations 



max I cot(it+a— v) cot(u— a— u)| = e 4V , min | cot(w+a— u) cot(u— a— u 



M=i 
or, 



U+e- 2H ± 1 



U+Te- 



-2H 



|u,|=l 

t/_e- 2H ± 1 



-2H 



,4V 



„4V 



(61) 



This curve is shown on Fig. [9] 

The ordered phases AIB2 and A2B1 are shown on Fig. [TQl 
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Figure 8: The boundary ordered and disordered regions in the (H, V)-plane for 
the homogeneous 6-vertex model with A = and U = cot(u) = 2 



The free energy per site in the disordered region is given by 

f(H, V) = \og{cos{u + a) cos(u — a)+ 

1 f Kl /cot(u + a)Lue- 2H + 1 cot(u - a)ue~ 2H + l\ duj 
2ni J_ Kl °^ \ uie~ 2H — cot(w + a) we~ 2H — cot(u — a) J uj 

1 r K ' 2 /cot(u + a)Loer 2H + 1 cot(u - a)ue- 2H + l\du 
2iriJ_ K2 ^ \ we~ 2H — cot(u + a) uie~ 2H — cot(u — a) / cj 



where if i 2 are defined by equations 

I cot(w + a — b) cot(u — a — b)\ — e 2V 
with cot(6) = e ±lK e- 2H . 
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Figure 9: The boundary between ordered and disordered regions in the (H, V)- 
plane for the inhomogeneous 6- vertex model with A = and U = cot (it) = 
3,T = cot(a) = 1 

8.3 Vertically non-homogeneous case 
8.3.1 

Assume the weights are homogeneous in the vertical direction and alternate in 
the horizontal direction with parameters alternating as ... , u + a, u — a, u + a, u — 
a, . . . . Positivity of weights in the region < u < ? requires < a < u. 
The eigenvalues of the transfer-matrix are given by the Bethe ansatz: 

A(u) = ((cos(« + a) cos(u - a)) N {-l) n e NH + {sm{u + a) sm{u - a)) N e - NH ) 

n 

e^ N - 2n ^ v l[cot(u-Vi) (63) 

i=l 

where Uj are distinct solutions to 

(cot(v - a) cot{v + a)) N = (-ly^er^ 11 
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Figure 10: Horizontally inhomogeneous ordered phases A\B 2 and A 2 Bi respec- 
tively 



The left side is the N-th power of 

cot(v) cot(a) + 1 — cot(w) cot(a) + 1 cot(v) 2 cot(a) 2 



1 



cot(a) — cot(v) — cot(a) — cot(u) cot(a) 2 — cot(v) 2 
From here it is easy to find the parametrization of solutions by roots of unity 

to cot(a) 2 er 2H + 1 



cot(w) z 



cot(a) 2 + oje~ 2H 



(64) 



where 



-iy 



8.3.2 Largest eigenvalue 

The analysis of the largest eigenvalue of the transfer-matrix in the limit N — > oo 
is similar to the previous cases. When 



max | cot (it — v)\ > e 

Ice? |= 1 



2V 



(65) 



where v and u> are related as in the absolute value of all factors cot(u — 
v)\e~ 2V is greater then one by the absolute value and the largest eigenvalue 
correspond to n = 0. The rest of the analysis is similar. 
Let us notations: 



T = cot(a), x = e 
Positivity of weights imply 



-H 



U = cot(u), 



cot(u) 



T > U > 

Proposition 1. The curve separating ordered phases from the disordered phase 
is 



X = 



YU+ - 1 
Y + U+ 



YU- - 
Y + U- 
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Proof. In the notations from above: 

, l2 .cot(it)cot(«) + l, 2 Us + lUs + 1 U 2 \s\ 2 + 2URe(s) + 1 
1 V n 1 cot(u) - cot(u) 1 8-U 8-U U 2 - 2URe(s) + \s\ 2 

The equation of the boundary of the disordered region is 

max | cot(u — v)\ = e 2V 



The equation defining v in terms of roots of unity lo can be solved explicitly 
for s — cot v: 

„2 



T 2 x 2 cu + 1 



T 2 + ujx 2 



Denote w = 



then 



2 _ (T 2 x 2 cos(a) + 1) + iT 2 x 2 sin(a) 
(T 2 + x 2 cos(a)) + ix 2 sin(a) 

, _ (T 2 x 2 cos(a) + l)(T 2 + x 2 cos(a)) + T 2 x 2 sin(a)x 2 sin(a) 



Re(s 2 ) 



Now, a simple algebra: 



H 4 = 



(T 2 + x 2 cos(a)) 2 + x 4 sin(a) 2 

T 4 x 4 + 2T 2 x 2 cos(a) + 1 
T 4 + 2T 2 x 2 cos(a) + x 4 



Re{s) 



Re{s 2 ) + \s\ 



As we vary uj along the unit circle, the function cot(it — v) has maximum at 
U)= 1. 

Taking this into account the equation for the boundary curve becomes 



C7V1 + x 2 T 2 ± VT 2 + x 2 



uVT 2 

Solving this equation for X — 

X = 



x 2 ± vT 



o2V 



2ff in terms of Y = e 2V we obtain 



YU- - 1 



Y + U- 



which is, after changing coordinates to X 1 ,Y 1 is the same curve as (|5T|) . This 
is one of the implications of the modular symmetry □ 

By modularity, i.e. by "rotating" the lattice with periodic boundary condi- 
tions by 90 degrees, all characteristics of the model with vertical inhomogeneities 
can be identifies with the corresponding characteristics of the model with hori- 
zontal inhomogeneities. 
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9 The free energy of the 6- vertex model 



The computation of the free energy for the 6- vertex model by taking the large 
TV, M limit of the partition function on a torus was outlines in section 17.11 In 
this section we will describe the free energy as a function of electric fields (H, V) 
and its basic properties. 

9.1 The phase diagram for A > 1 
9.1.1 

The weights a, b, and c in this region satisfy one of the two inequalities, either 
a > b + c or b > a + c. 

If a > b + c, the Boltzmann weights a, 6, and c can be parameterized as 

a = r sinh( A + r)), b = r sinh(A) ,c = r sinh(r?) (66) 

with A, r\ > 0. 

If a + c < b, the Boltzmann weights can be parameterized as 

a = r sinh( A — r)),b = r sinh(A) ,c = r sinh(r?) (67) 

with < rj < A. 

For both of these parametrization of weights A = cosh(^) . 

The phase diagram of the model for a > b + c (and, therefore, a > b) is 
shown on Fig. [T21and for b > a + c (and, therefore, a < b) on Fig. IT51 

9.1.2 

When magnetic fields (if, V) are in one of the regions Ai, Bi of the phase dia- 
gram, the system in the thermodynamic limit is described by the translationally 
invariant Gibbs measure supported on the corresponding frozen (ordered) con- 
figurations. There are four frozen configurations A±, A%, B%, and B2, shown on 
Fig. [TTJ For a finite but large grid the probability of any other state is of the 
order at most exp(— aN) for some positive a. 

Local correlation functions in a frozen state are products of expectation 
values of characteristic functions of edges. 

lim (cr ei . . . er e „)jv = Cei (S) . . ■ a Cn (S) 
where 5 is the one of the ferromagnetic states Ai , Bi . 
9.1.3 

The boundary between ordered phases in the (H, V^-plane and disordered phases, 
as in the free fermionic case, is determined by the next to the largest eigenvalue 
of the row-to-row transfer matrix. 
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Without going into the details of the computations as we did in the free 
fermionic case we will just give present answers. 

• a > b + c, see Fig. Hlfl 

Ai-region: V + H > 0, cosh(2ff) < A, 

{e 2H ~b/a)(e 2V -b/a) > (c/a) 2 , e 2H > b/a, cosh(2iJ) > A, 
^2-region: V + H < 0, cosh(2i7) < A, 

(e~ 2H -b/a)(e- 2V -b/a) > (c/a) 2 , e~ 2ff > b/a, cosh(2H) > A, 
Si-region: (e 2H - a/b){e- 2V - a/b) > (c/b) 2 , e 2H > a/b, 
B 2 -region: {e' 2 " - a/b)(e 2V - a/b) > (c/b) 2 , e' 211 > a/b. 

• b > a + c, see Fig. H31 

^-region: {e 2H - b/a)(e 2V - b/a) > (c/a) 2 , e 2H > b/a; 
^2-region: (e" 2H - b/a){e~ 2V - b/a) > (c/a) 2 , e~ 2H > b/a; 
Bi-region: 7 - H > 0, cosh(2ff) < A, 

(e 2H -a/b)(e~ 2V -a/b)>(c/b) 2 , e 2H > a/b, cosh(2if) > A, 
B 2 -region: V - H < 0, cosh(2i7) < A, 

(e' 2H -a/6)(e 21/ -a/6)> (c/6) 2 , e - 2H > a/6, cosh(2if) > A. 

The free energy is a linear function in H and V in the four frozen regions: 

f = -1na-H-V in A 1 , 

f = -\nb + H -V in B 2) (68) 

/ = -hxa + H + V in A%, 

f = -l n b-H + V in Si. 

The regions D\ and D2 are disordered phases. If (H, V) is in one of these 
regions, local correlation functions are determined by the unique Gibbs measure 
with the polarization given by the gradient of the free energy. In this phase the 
system is disordered, which means that local correlation functions decay as a 
power of the distance d(a,ej) between a and ej when d(ei,ej) — > 00. 



Figure 11: Four frozen configurations of the ferromagnetic phase 



4G 




Figure 12: The phase diagram in the (H, F)-plane for a = 2, b = 1, and c = 0.8 

In the regions D\ and D2 the free energy is given by |37j : 

/(if, 7) - min(min ( E 1 - H - (1 - 2a) V - -L / ln(- - , 
q \ 27TI a ao — a z z / 

min (£ 2 + if - (1 - 2a)V - -L / m(^— — + - - „ (69) 
" \ 27TJ ao ao — a z J 

where p[z) can be found from the integral equation 

/ \ 1 1 f P( w ) , 1 f P( w ) , 

= - + tt- r dtfl " ^ / r dto, (70) 

in which 

^i(w) = — , z 2 {w) = + 2A. 

2A — w w 

p[z) satisfies the following normalization condition: 

1 f 

a = / p(z)dz. 

The contour of integration C (in the complex z-plane) is symmetric with respect 
to the conjugation z — > z, is dependent on H and is defined by the condition 
that the form p(z)dz has purely imaginary values on the vectors tangent to C: 



Re(p(z)dz) 



= 0. 

26 c 
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Figure 13: The phase diagram in the (H, F)-plane for a = 1, b = 2, and c = 0.8 



The formula (|69p for the free energy follows from the Bethe Ansatz diagonal- 
ization of the row-to-row transfer-matrix. It relies on a number of conjectures 
that are supported by numerical and analytical evidence and in physics are 
taken for granted. However, there is no rigorous proof. 

There are two points where three phases coexist (two frozen and one dis- 
ordered phase). These points are called tricritical. The angle 6 between the 
boundaries of D\ (or D2) at a tricritical point is given by 



v ' c 2 + 2min(a,&) 2 (A 2 - 1)' 

The existence of such points makes the 6- vertex model (and its degeneration 
known as the 5-vertex model |15j ) remarkably different from dimer models [20j 
where generic singularities in the phase diagram are cusps. Physically, the 
existence of singular points where two curves meet at the finite angle manifests 
the presence of interaction in the 6- vertex model. 

Notice that when A = 1 the phase diagram of the model has a cusp at the 
point H = V = 0. This is the transitional point between the region A > 1 and 
the region |A| < 1 which is described below. 
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Figure 14: The phase diagram in the (H, T^)-plane for a = 1, b = 2, and c = 2 
9.2 The phase diagram |A| < 1 

In this case, the Boltzmann weights have a convenient parametrization by 
trigonometric functions. When 1 > A > 1 

a = r sin(A — 7), b = r sin(A), c — r sin (7), 

where < 7 < n/2, 7 < A < 7r, and A = cos 7. 
When > A > -1 

a = r sin(7 — A), b = r sin(A), c = r sin(7), 

where < 7 < w/2, tt — 7 < A < tt, and A = — cos 7. 

The phase diagram of the 6-vertex model with |A| < 1 is shown on Fig. 
[TH The phases Ai , Bi are frozen and identical to the frozen phases for A > 1 . 
The phase D is disordered. For magnetic fields (H, V) the Gibbs measure is 
translationally invariant with the slope (h,v) = ( ^S^ , d ^Qy V * > )■ 
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The frozen phases can be described by the following inequalities: 

Ai-region: (e 2H - b/a)(e 2V - b/a) > (c/a) 2 , e 2H > b/a, 

A 2 -region: (e- 2H - b/a)( e - 2V - b/a) > (c/a) 2 , e~ 2H > b/a, 

Bx-region: (e 2H - a/b)( e - 2V - a/b) > (c/b) 2 , e 2H > a/b, (71) 

5 2 -region: (eT 2H - a/b)(e 2V - a/b) > (c/b) 2 , e' 2 " > a/b. 

The free energy function in the frozen regions is still given by the formulae 
(|69[) . The first derivatives of the free energy are continuous at the boundary of 
frozen phases, The second derivative is continuous in the tangent direction at 
the boundary of frozen phases and is singular in the normal direction. 

It is smooth in the disordered region where it is given by (I69[) which, as in 
case A > 1 involves a solution to the integral equation (TfOl) . The contour of 
integration in ([70|) is closed for zero magnetic fields and, therefore, the equation 
(|70[) can be solved explicitly by the Fourier transformation [3] . 

The 6- vertex Gibbs measure with zero magnetic fields converges in the ther- 
modynamic limit to the superposition of translationally invariant Gibbs mea- 
sures with the slope (1/2, 1/2). There are two such measures. They correspond 
to the double degeneracy of the largest eigenvalue of the row-to-row transfer- 
matrix [3]. 

There is a very interesting relationship between the 6-vertex model in zero 
magnetic fields and the highest weight representation theory of the correspond- 
ing quantum affine algebra. The double degeneracy of the Gibbs measure with 
the slope (1/2, 1/2) corresponds to the fact that there are two integrable irre- 
ducible representations of sZ 2 at level one. Correlation functions in this case can 
be computed using q- vertex operators [T7] . For latest developments see [7] . 



9.3 The phase diagram A < — 1 
9.3.1 The phase diagram 

The Boltzmann weights for these values of A can be conveniently parameterized 
as 

a = r sinh(f7 — A), b = r sinh(A), c = r sinh^), (72) 

where < A < i] and A = — cosh 77. 

The Gibbs measure in thermodynamic limit depends on the value of mag- 
netic fields. The phase diagram in this case is shown on Fig. [l5]for b/a > 1. In 
the parameterization (IT2"j) this correspond toO<A<?7/2. When 7//2 < A < 7/ 
the 4-tentacled "amoeba" is tilted in the opposite direction as on Fig. [12J 

When (H, V) is in one of the A i7 Bi regions in the phase diagram the Gibbs 
measure is supported on the corresponding frozen configuration, see Fig. 1111 

The boundary between ordered phases Ai , Bi and the disordered phase D is 
given by inequalities (|7ip . The free energy in these regions is linear in electric 
fields and is given by (|69|) . 

If (H, V) is in the region D, the Gibbs measure is the translationally invariant 
measure with the polarization (ft,, v) determined by (|54[) . The free energy in this 
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Figure 15: The phase diagram in the (H, T^)-plane for a = 1, b = 2, and c = 6 

case is determined by the solution to the linear integral equation ([70)1 and is 
given by the formula 

If (H, V) is in the region A, the Gibbs measure is the superposition of two 
Gibbs measures with the polyarization (1/2, 1/2). In the limit A — > — oo these 
two measures degenerate to two measures supported on configurations C\, Ca, 
respectively, shown on Fig. [TH1 For a finite A the support of these measures 
consists of configurations which differ from C\ and C*2 in finitely many places 
on the lattice. 

Remark 4. Any two configurations lying in the support of each of these Gibbs 
measures can be obtained from C\ or C2 via flipping the path at a vertex "up" 
or "down" as it is shown on Fig. finitely many times. It is also clear that it 
takes infinitely many flips to go from C\ to Ci. 

9.3.2 The antiferromagnetic region 

The 6-vertex model in the phase A is disordered and is also noncritical. Here 
the non-criticality means that the local correlation function (a ei a ej ) decays as 
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exp(— ad(ei, ej)) with some positive a as the distance d(ei,ej) between and 
ej increases to infinity. 

The free energy in the A-region can be explicitly computed by solving the 
equation (|T0|) . In this case the largest eigenvalue will correspond to n = N/2 , 
the contour of integration in (|70p is closed, and the equation can be solved by 
the Fourier transform. 0. 

The boundary between the antiferromagnetic region A and the disordered 
region D can be derived similarly to the boundaries of the ferromagnetic regions 
Ai and Bi by analyzing next to the largest eigenvalue of the row-to-row transfer- 
matrix. This computation was done in [37], [21]. The result is a simple closed 
curve, which can be described parameterically as 



where 



H(s) = E(s), V(s) = E(r) 
E(ip) = cosh - 



1 



and 



dn{ty\l-v) 
\s\<2rj, 

1 + max(6/a, a/b)e v 



max(6/a, a/b) + e n 
The parameter v is defined by the equation r\K{y) = ttK'(p), where 

fir/2 

K(v)= (l-^sin 2 (60r 1/2 d<9 K'{y) = / (1- (1- v) sm 2 {9))-^ 2 d6. 
Jo Jo 

The curve is invariant with respect to the reflections (H, V) — > (-H, —V) 
and (H, V) — > (V, H) since the function S satisfies the identities 

3(y>) = -E(-<p), E(r]-<p) =E(r) + <p). 



"Strictly speaking this is a conjecture supported by the numerical evidence 



it 





































n r 





Figure 16: The configurations C\ and C2. 
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-J*-> 

^ down 

-FF^ -F- 

Figure 17: The elementary up and down fluctuations in the antiferromagnetic 
phase. 

This function is also 4?7-periodic: S(4ry + ip) = S(y>). 

As it was shown in |29j this curve is algebraic in e H and e v and can be 
written as 

({l-v cosh 2 Vo) cosh 2 H + sinh 2 Vq-(1-u) cosh Vb cosh H cosh = 

(l-i/ cosh 2 Vo) sinh 2 Vb cosh 2 V sinh 2 H{l-v cosh 2 if), (73) 

where Vb is the positive value of 1/ on the curve when H = 0. Notice that i/ 
depends on the Boltzmann weights a, 6, c only through r\. 



10 Some asymptotics of the free energy 

10.1 The scaling near the boundary of the D-region 

Assume that Hq = (Hq,Vq) is a regular point at the boundary between the 
disordered region and and the A\ -region, see FigfT4l Recall that this boundary 
is the curve defined by the equation 

g(H, V)=0 

where 

g(H, V) = Hb/a + e2 Pf b/a ) - 2V. (74) 

Denote the normal vector to the boundary of the Z?-region at Hq by n and 
the tangent vector pointing inside of the region D by r. 

We will study the asymptotic of the free energy along the curves H(r, s, t) — 
Ho + r 2 sn + rtr, as r — > 0. It is clear that H (r, s, t) is in the D-region if s > . 

Theorem 2. Let H(r,s,t) be defined as above. The asymptotic of the free 
energy of the 6-vertex model in the limit r — > is given by 

f(H(r, s, t)) - f lm (H(r, s, t)) + v (s, t)r 3 + 0(r 5 ), (75) 
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where fu n (H, V) = — ln(a) — H — V and 

r)(s,t)=-K(6s + t 2 ) 3/2 . (76) 

Here the constants k and 9 depend on the Boltzmann weights of the model and 
on (Hq, Vq) and are given by 

K=^d 2 H g(H ,Vo) 

and 

4+(d H g(H ,V )) 2 
2d%g(H ,V ) ' 

where g(H, V) is defined in J7^[ ). 

Moreover, djjg(Ho,Vo) > and, therefore, 9 > 0. 

We refer the reader to [29] for the details. This behavior is universal in a 
sense that the exponent 3/2 is the same for all points at the boundary. 

10.2 The scaling in the tentacle 

Here, to be specific we assume that a > b. The theorem below describes the 
asymptotic of the free energy function when H — > +oo and 

\ ^ ~ &b e ~ 2H ^ A ^ + &b e ~ 2H > H ^°°- ^ 
These values of (H, V) describe points inside the right "tentacle" on the Fig. 

Let us parameterize these values of V as 

^=iln(6/a) + ^e-- 

where f3 £ [-1, 1]. 

Theorem 3. When H — > oo and j3 € [—1,1] the asymptotic of the free energy 
is given by the following formula: 

f(H,V) = -hn(ab)-H- 

^re- 2H (f3+ -yi^--/3arccos(/3)) +0(e- 4H ), 

lab V 7T 7T / 

The proof is given in [29] . 
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10.3 The 5-vertex limit 

The 5-vertex model can be obtained as the limit of the 6-vertex model when 
A — > oo. Magnetic fields in this limit behave as follows: 

• a > b + c. In the parametrization (1661) after changing variables H = ^ + 1, 
and V = — \ + m take the limit 77 — > 00 keeping A fixed. The weights will 
converge (up to a common factor) to: 

ai : a 2 : h : b 2 : a : c 2 -> e X+l+m : e*" 1 -™ : ( e A - e - A ) e <-™ : : 1 : 1 

• a + c < b. In the parametrization (|67"1) after changing variables H = 3 + 1, 
and V = T; + m take the limit 77 — > 00 keeping £ = A — 77 fixed. The weights 
will converge (up to a common factor) to: 

a x : a 2 : h : b 2 : c x : c 2 -> (e« - e^)e' +m : : e i+l ~ m : e^' +m : 1 : 1 

The two limits are related by inverting horizontal arrows. From now on we will 
focus on the 5-vertex model obtained by the limit from the 6-vertex one when 
a > b + c. 

The phase diagram of the 5-vertex model is easier then the one for the 
6-vertex model but still sufficiently interesting. Perhaps the most interesting 
feature is that the existence of the tricritical point in the phase diagram. 

We will use the parameter 

7 = e- 2X 

Notice that 7 < 1. 

The frozen regions on the phase diagram of the 5-vertex model, denoted on 
Fig. [18] as Ai, A 2 , and B\, can be described by the following inequalities: 

Ai-region: m > —I, I < 0, 

e 2m > 1-7(1 -e- 21 ), 1>1; 
A 2 -region: m < -I, I < 0, (78) 

e 2m <l--(l-e- 21 ), 1>1; 
7 

1 w_-2m 1 s ^ 7 ox ^ 1 



^-region: {e M - ){e~ Zm - ) > - ' „ , e M > 



1— 7 1 — 7 (1 — 7) 2 ' 1— 7' 

As it follows from results [15] the limit from the 6-vertex model to the 5- 
vertex model commutes with the thermodynamical limit and for the free energy 
of the 5-vertex model we can use the formula 

/ 5 (Z,m)= lim (/(r ? /2 + /,-7 7 /2 + m)-/(7 ? /2 ) -r ? /2)), (79) 
where f(H, V) is the free energy of the 6-vertex model. 
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Figure 18: The phase diagram of the 5-vertex model with 7 = 1/4 (/? = e 2h ). 

10.4 The asymptotic of the free energy near the tricritical 
point in the 5-vertex model 

The disordered region D near the tricritical point forms a corner 

— l + 0(l 2 ) < m< --yl + 0(l 2 ), h^0 + . 

7 

The angle 9 between the boundaries of the disordered region at this point is 
given by 

cos(0) = 

1 + 7^ 

One can argue that the finiteness of the angle 9 manifests the presence of 
interaction in the model. In comparison, translation invariant dimer models 
most likely can only have cusps as such singularities. 

Let 7 < k < - and 

' — — 7 

m = —kl, 

As it was shown in [6 for m — —I and in 29 for m = —kl with 7 < k < - the 
asymptotic of the free energy as I +0 is given by 

f(l,-kl) = Cl (k,j)l + c 2 (k, 1 )l 5 / 3 + 0{l 7 '*). (80) 

where 

c 1 (fc,7) = T -^(-(l + fc)(l+7)+4^), (81) 
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and _ 

, (l|T) = w ,«l!W^, (82) 

The scaling along any ray inside the corner near the tricritical point in the 
6-vertex model differ from this only by in coefficients. The exponent ft 5 / 3 is the 
same. 

10.5 The limit A ->■ -l - 

If A = — 1, the region A consists of one point located at the origin. 

It is easy to find the asymptotic of the function ^(ip) when r\ — > 0+, or 
A — > — 1 — . In this limit K' — > 7r/2, K — > Using the asymptotic dn ( M |\_ m ) ~ 
1 + i(l — m)sin 2 (u) when m — > 1—, and cosh _1 (a;) <~ ±^2(x — 1), when 



2 ( 

1+, we have 



and 



\b-a\r] 
i o = —r— ■ 



E(<p) ~ 4e 2 i sin( — (^). 
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The gap in the spectrum of elementary excitations vanishes in this limit at 

the same rate as 2(93). At the lattice distances of order e" 5 ^ the theory has a 
scaling limit and become the relativistic SU(2) chiral Thirring model. Corre- 
lation function of vertical and horizontal edges in the 6-vertex model become 
correlation functions of the currents in the Thirring model. 

10.6 The convexity of the free energy 

The following identity holds in the region D [4].|26j: 

fH.nfv.v - Ih,v = f — ) • (^3) 

Here g = 75771. The constant Dq does not vanish in the D- region including 
its boundary. It is determined by the solution to the integral equation for the 
density p{z). 

Directly from the definition of the free energy we have 
f HH = hm <^)-n { R)f> 

■ N.M^oo NM 

where n(L) and n(R) are the number of arrows pointing to the left and the 
number of arrows pointing to the right, respectively. 

Therefore, the matrix didjf of second derivatives with respect to H and V 
is positive definite. 
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As it follows from the asymptotical behavior of the free energy near the 
boundary of the D-phase, despite the fact that the Hessian is nonzero and finite 
at the boundary of the interface, the second derivative of the free energy in the 
transversal direction at a generic point of the interface develops a singularity. 

11 The Legendre Transform of the Free Energy 

The Legendre transform of the free energy 

sup(xH + yV + f(H, V)) 

as a function of (x, y) is defined for — 1 < x, y < 1. 

The variables x and y are known as polarizations and are related to the slope 
of the Gibbs measure as x = 2h — 1 and y = 2v — 1. We will write the Legendre 
transform of the free energy as a function of (ft, v) 

a(h,v) = sup((2h-l)H+(2v-l)V + f{H,V)). (84) 

H,V ^ ' 

a(h, v) is defined on < ft, v < 1. 

For the periodic boundary conditions the surface tension function has the 
following symmetries: 

a(x, y) = a(y, x) = a(-x, -y) = a(-y, -x). 

The last two equalities follow from the fact that if all arrows are reversed, a is the 
same, but the signs of x and y are changed. It follows that 07, (ft, v) — a v (v, h) 
and a v (h,v) = ah(v,h). 

The function f(H, V) is linear in the domains that correspond to conic and 
corner singularities of a. Outside of these domains (in the disordered domain 
D) we have 

Vcr o V/ = id D , V/o V<7 = id V /(£>). (85) 

Here the gradient of a function as a mapping K 2 — > M 2 . 

When the 6-vertex model is formulated in terms of the height function, the 
Legendre transform of the free energy can be regarded as a surface tension. The 
surface in this terminology is the graph of of the height function. 

11.1 

Now let us describe some analytical properties of the function <r(ft, v) is obtained 
as the Legendre transform of the free energy. The Legendre transform maps the 
regions where the free energy is linear with the slope (±1, ±1) to the corners of 
the unit square D = {(ft, v)\ < ft < 1, < v < 1}. For example, the region 
A\ is mapped to the corner h = 1 and v = 1 and the region Bi is mapped to 
the corner h = 1 and v = 0. The Legendre transform maps the tentacles of the 
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disordered region to the regions adjacent to the boundary of the unit square. 
For example, the tentacle between A\ and B\ frozen regions is mapped into a 
neighborhood of ft = 1 boundary of D, i.e. ft, 1 and < v < 1. 

Applying the Legendre transform to asymptotics of the free energy in the 
tentacle between A\ and B\ frozen regions we get 

H(h,v) = -im (Z^ 5 J^_) , v M = \w/a) + l (I-Mcotwi- 

and 

(7(fe,t;) = (l-/i)lnf^ . /T^ ) -(l-h) + vHb/a)-\n(b), (86) 
\ c z sm(7r(l — v)j J 

Here ft, 1— and < v < 1. From ([55]) we see that a(l,v) = vhx(b/a) — hi(6), 
i.e. cr is linear on the boundary h = 1 of D. Therefore, its asymptotics near the 
boundary ft, = 1 is given by 

a(h, v)=v hx(b/a) - hx(b) + (1 - ft) ln(l - ft) + 0(1 - ft), 

as ft —> 1— and < v < 1. We note that this expansion is valid when (1 — 
ft)/sin(7r(l - u)) < 1. 

Similarly, considering other tentacles of the region D, we conclude that the 
surface tension function is linear on the boundary of D. 



11.2 

Next let us find the asymptotics of a at the corners of D in the case when all 
points of the interfaces between frozen and disordered regions are regular, i.e. 
when A < 1. We use the asymptotics of the free energy near the interface 
between A\ and D regions (|75|) . 

First let us fix the point (Hq, Vq) on the interface and the scaling factor r in 
(|75|) . Then from the Legendre transform we get 

4 {d H g) 2 +4 
l-ft 6d H g + ^rt 



and 



It follows that 



l-v 2{-6 + rtd B g) 
In the vicinity of the boundary r and, hence, 

l-ft _ _d H g _ 1 - ft/a e~ 2V ° 
l-v ~ 2~ _ 1 - b/a e~ 2H ° 



(87) 
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as h, v — > 1. Thus, under the Legendre transform, the slope of the line which 
approaches the corner h — v — 1 depends on the boundary point on the interface 
between the frozen and disordered regions. 

It follows that the first terms of the asymptotics of a at the corner h = v = 1 
are given by 

a(h,v) = -lna - 2(1 - h)H (h,v) - 2(1 - v)V (h,v), 

where H (h, v) and Vo(h, v) can be found from (|87)) and g(H 0} Vo) = 0. 

When | A | < 1 the function a is strictly convex and smooth for all < h, v < 
1. It develops conical singularities near the boundary. 

When A < —1, in addition to the singularities on the boundary, a has a 
conical singularity at the point (1/2,1/2). It corresponds to the "central flat 
part" of the free energy /, see Fig. [15] 

When A > 1 the function a has corner singularities along the boundary as 
in the other cases. In addition to this, it has a corner singularity along the 
diagonal v — h if a > b and v = 1 — h ii a < b. We refer the reader to [5] for 
further details on singularities of a in the case when A > 1. 

When A = — 1 function a has a corner singularity at ((1/2, 1/2). 

12 The limit shape phenomenon 

12.1 The Height Function for the 6-vertex model 

Consider the square grid L c C M 2 with the step e. Let D C R 2 be a domain 
in R 2 . Denote by D e a domain in the square lattice which corresponds to the 
intersection DO L e , assuming that the intersection is generic, i.e. the boundary 
of D does not intersect vertices of L e . 

Faces of D e which do not intersect the boundary 3D of D are called inner 
faces. Faces of D e which intersect the boundary of D are called boundary faces. 

A height function h is an integer-valued function on the faces of D e of the 
grid Lm (including the outer faces) which is 

• non-decreasing when going up or to the right, 

• if fi and fi are neighboring faces then /i(/i) — h(f2) = —1, 0, 1. 

Boundary value of the height function is its restriction to the outer faces. 
Given a function on the set of boundary faces denote ^H(h^) the space of 
all height functions with the boundary value h^°\ Choose a marked face /o at 
the boundary. A height function is normalized at this face if h(fo) = 0. 

Proposition 2. There is a bijection between the states of the 6-vertex model 
with fixed boundary conditions, and height functions with corresponding bound- 
ary values normalized at fo . 
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Figure 19: The values of the height function for the configuration of paths given 
on Fig. [5] 

Proof. Indeed, given a height function consider its "level curves," i.e. paths 
on D e , where the height function changes its value by 1, see Fig. Q1JJ Clearly 
this defines a state for the 6- vertex model on D e with the boundary conditions 
determined by the boundary values of the height function. 

On the other hand, given a state in the 6-vertex model, consider the corre- 
sponding configuration of paths. That there is a unique height function whose 
level curves are these paths and which satisfies the condition h = at /o 

It is clear that this correspondence is a bijection. □ 

There is a natural partial order on the set of height functions with given 
boundary values. One function is bigger than the other, hi < hi, if it is entirely 
above the other, i.e. if h\(x) < hi(x) for all x in the domain. There exist 
the minimum ft. m i n and the maximum /i max height functions such that h mm < 
h < /imax for all height functions h. for DW boundary conditions maximal and 
minimal height functions ar shown on Fig. 1201 

The characteristic function of an edge is related to the hight function as 

a e = h(f+) - h(f-) (88) 

where faces ff- are adjacent to e, /+ is to the right of e for the vertical edge 
and on the top for the horizontal one. The face /~ is on the left for the vertical 
edge and below e for the horizontal edge. 

Thus, we can consider the 6-vertex model on a domain as a theory of fluctu- 
ating discrete surfaces constrained between h m - m and ft. max - Each surface occurs 
with probability given by the Boltzmann weights of the 6-vertex model. 

The height function does not exist when the region i snot simply-connected 
or when a square lattice is on a surface with non-trivial fundamental group. One 
can draw analogy between states of the 6-vertex model and 1-forms. States on a 
domain with trivial fundamental group can be regarded as exact forms uj = dh 
where h is a height function. 
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12.2 Stabilizing Fixed Boundary Conditions 

Recall that the height function is a monotonic integer-valued function on the 
faces of the grid, which satisfies the Lipschitz condition (it changes at most by 
1 on any two adjacent faces). 

The normalized height function is a piecewise constant function on D e with 
the value 

h noim {x iy ) = eh e (n,m). 

where h e (n, m) is a height function on D e . 

Normalized height functions on the boundary of D satisfy the inequality 

\h(x,y)-h(x',y')\ < \x-x'\ + \y-y'\. 

As for non-normalized height functions there is a natural partial ordering on 
the set of all normalized height functions with given boundary values: hi > h-2 
if hi(x) > h2(x) for all x E D. We define the operations 

hi V /12 = min(/ii(x), h2(x)), hi A /12 = max(/ii(x), h2{x)). 
It is clear that 

h\ V h<i < hi , h-z < hi A /i2 

It is also clear that in this partial ordering there is a unique minimal and max- 
imal height functions, which we denote by h mm and h max , respectively. 

The boundary value of the height function defines a piecewise constant func- 
tion on boundary faces of D t . Consider a sequence of domains D t with e — > 0. 
Stabilizing fixed boundary conditions for height functions is a sequence {h^} 
of functions on boundary faces of D e such that h^ — > M°) where h^ is a 
continuous function of dD. 

It is clear that in the limit e — > we have at least two characteristic scales. 
At the macro scopical scale we can " see" the region Del 2 . Normalized height 
functions in the limit e — > will be functions on D. As we will see in the next 




Figure 20: The minimum and maximum height functions for the DW boundary 
conditions. 
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section, the height function in the 6- vertex model develops deterministic limit 
shape ho(x,y) at this scale. It is reflected in the structure of local correlation 
functions. Let a be edges with coordinates (n^m^) = (— , When e — > 
and (xj, yi) are fixed 



The randomness remain at the smaller scale, and in particular at the lattice, 
microscopical, scale. If coordinates of edges are {n^mi) — (|+Arij, |+Amj) 
we expect that in the limit e — > correlation functions have the limit 



where the correlation function at the right side is taken with respect to the trans- 
lation invariant Gibbs nmeasure with the average polarization (d x ho(x, y), d y ho(x, y)). 
The correlator in the r.h.s. depends on the polarization and on An^— Arij , Ami — 
Arrij, i.e. it is translation invariant. 

12.3 The Variational Principle 
12.3.1 

Here we will outline the derivation of the variational problem which determines 
the limit shape of the height function. 

First, consider the sequence of rectangular domain of size M, M when N, M — > 
oo such that a — N/M is finite and the boundary values of the height function 
stabilize to the function (f>(x,y) — hy + vx,x £ [0, a],y £ [0,1]. The partition 
function of such system has the asymptotic 



where a(h, v) is the Legendre transform of the free energy for the torus. 

For a domain D e , choose a subdivision of it into a collection of small rect- 
angles. Taking into account ([89]) . the partition function of the 6- vertex model 
with zero electric fields and stabilizing boundary conditions can be written as 



h 

When e — > the size of \D e \ the region is increasing and the the leading contribu- 
tion to the sum comes from the height function which minimizes the functional 



One can introduce the extra weight q vol ( h ) to the partition function of the 
6- vertex model. It corresponds to inhomogeneous electric fields [29]. If e — > 
and q = exp(— Ae) the leading contribution to the partition function comes from 
the height function which minimizes the functional 



< er ei . . .<j ek >= d il ho{xx,yi) ■ ■■d ik ho(x k ,yk) + 0(e <=) 



< a ei . . .a Ek >-> a ei ... a &k >(d x h (x, v ),d y h (x, v )) 



Zn,m oc exp(NMa(h,v)) 



(89) 
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12.3.2 The variational principle 

Thus, in order to find the limit shape in the thermodynamical limit of the 
6- vertex model on a disc with Dirichlet boundary conditions (fo, we should 
minimize the functional 

I x [<p] = I cr(\7ip)d 2 x + X I ipd 2 x, (90) 

J D J D 

on the space L(D,(po) of functions satisfying the Lipshitz condition 

\ip(x,y) - f(x',y')\ <\x- x'\ + \y - y'\ 
and the boundary conditions monotonically increasing in x and y directions 

<f\dD = ¥>o- 

Since a is convex, the minimizer is unique when it exists. Thus, we should 
expect that the variational problem (|90[) has a unique solution. 

The large deviation principle applied to this situation should result in the 
convergence in probability, as e — > of random normalized height functions 
h(x,y) to the minimizer of (|90[) . 

12.3.3 

If the vector \7h(x,y) is not a singular point of tr, the minimizer h satisfies the 
Euler-Lagrange equation in a neighborhood of (x,y) 

div(Ver o Vh) = A. (91) 

We can also rewrite this equation in the form 

Va(Vh(x, y)) = ^ (x, y) + (-g y (x, y), g x (x, y)), (92) 

where g is an unknown function such that g xy (x, y) = g yx (x, y). It is determined 
by the boundary conditions for h. 
Applying (pj|) . it follows that 

Vh(x,y) = Vf(-x-g y [x,y),^y + g x (x,yf). (93) 

From the definition of the slope, see (|54l) . it follows that \Jh\ < 1 and 
|/v| < 1- Thus, if the minimizer h is differentiable at (x,y), it satisfies the 
constrains \h x \ < 1 and \h y \ < 1. It is given that h x ,h y > 0, hence, < h x < 1 
and < h y < 1. 

In particular, one can choose g(x, y) = 0. In this case the function is the 
minimizer of the rate functional /a[^] with very special boundary conditions. 
For infinite region D (and finite A) this height function reproduces the free 
energy as a function of electric fields. 
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The limit shape height function is a real analytic function almost everywhere. 
One of the corollaries of (|93[) is that curves on which it is not analytic ( interfaces 
of the limit shape) are images of boundaries between different phases in the 
phase diagram of f(H, V). The ferroelectric and antiferroelectric regions, where 
/ is linear, correspond to the regions where h is linear. 

From (1931) and the asymptotic of free energy near boundaries between dif- 
ferent phases, we can make some general conclusions about the the structure of 
limit shapes. 

First conclusion is that near regular pieces of boundary between regions 
where the height function is real analytic it behaves as h(x, y) — h(xo,yo) — d 3 ^ 2 
where d is the distance from (x, y) to the closest point (xq, yo) at the boundary. 

Second obvious conclusion is that if (x,y) is inside the corner singularity 
in the boundary between smooth parts, the height function behaves near the 
corner as h(x,y) — hrXo,yo) — d 5 ^ 3 where d is the distance from (x, y) to the 
corner, and (xo,yo) are coordinates of the corner. 

12.4 Limit shapes for inhomogeneous models 

So far the analysis of limit shapes was done in the homogeneous 6-vertex model. 
The variational problem determining limit shapes involves the free energy as the 
function of the polarization. The extension of the analysis outlined above to an 
inhomogeneous periodically weighted case is straightforward. The variational 
principle is the same same but the function a is determined from the Bethe 
ansatz for the inhomogeneous model. It is again the Legendre transform of the 
free energy as a function of electric fields. 

The computation of the free energy from the large N asymptotic of Bethe 
ansatz equations is similar and based on similar conjecture about the accumu- 
lation of solutions to Bethe equations on a curve. In the free fermionic case it 
is illustrated in section [8l If the size of the fundamental domain is k x m one 
should expect 2{k + m) cusps in the boundary between ordered and disordered 
regions. 

12.5 Higher spin 6-vertex model 

The weights in the 6-vertex model can be identified with matrix elements of the 
i?-matrix of U q (sl2) in the tensor product of two 2-dimensional representations 
of this algebra. We denoted such matrix Ry-^>{z). 

In a similar fashion one can consider matrix elements of Ry 1,l2 \z) as weights 
of the higher spin generalization of the 6-vertex model. We will call it higher 
spin 6-vertex model. 

A state in such models is an assignment of a weight of the irreducible rep- 
resentation V^ 1 ^ to every horizontal edge, and of a weight of the irreducible 
representation of to every vertical edge. 

Natural local observables in such model are sums of products of spin func- 
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tions of an edge: 

s e (S) = the weight s assigned to e in the state S, 
Spin variables define the height function locally as 

s e = h(f+) - h(f-) 

where ff~ are faces adjacent to e. When the domain has trivial fundamental 
group the local height function extends to a global one. Otherwise one should 
the region into simply-connected pieces. 
The height function has the property 

\h(n, m) — h(n' , m')| < li\n — n'\ + Ja|m — m\, 

Given a sequence of domains D e with e — > one should expect that the 
normalized height function eh(en,em) develops the limit shape ho(x,y) which 
minimizes the functional 

I^M[h] = [ a^ h ' h \Vh)d 2 x 

J D 

subset to the Dirichlet boundary conditions and the constraint 
\h(x, y) - h(x', 2/01 < ills - A + h\v - y'\, 

The properties of such model and of its free energy as function of electric 
fields is an interesting problem which need further research. 

13 Semiclassical limits 

13.1 The semiclassical limits in Bethe states 

The spectrum of quantum spin chains described in section U in terms of Bethe 
equations has a natural limit when rrii = Ri/h, q — e , and h —> with fixed 
Ri. In this limit the quantum spin chain with the quantum monodromy matrix 
(135[) becomes the classical spin chain with the monodromy matrix 

Semiclassical eigenvectors of quantum transfer-matrices correspond to solu- 
tions to Bethe equations which accumulate along contours representing branch 
cuts on the spectral curve of T(z), see [31] and for more recent results [36] . 

13.2 The semiclassical limits in the higher spin 6-vertex 
model 

Relatively little known about the semiclassical limit of the higher spin 6-vertex 
model. This is the limit when l\ = Ri/h, and h — > 0. Here, depending on the 
values of A, h = r\ or h = 7. 
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The first problem is to find the asymptotic of the conjugation action of the 
i?-matrix in this limit. The mapping 

x R(luh)( u <) xR (h,l2)( u yi 

is an automorphism of En^V 1 ^ 1 ^ C3> V^) which as h — > becomes the Poisson 
automorphism p( Rl,R2 \u) of the Poisson algebra of functions on S^ Rl ' €5 S^ R2 \ 

When A = 1 this automorphism was computed in [35] . It is easy to extend 
this results for A ^ 1. For constant i?-matrices see [3D]. 

Next problem is to find the semiclassical asymptotic for the i?-matrix con- 
sidered as an "evolution operators". For this one should choose two Lagrangian 
submanifolds in S^ Rl ^ ® S^ R2 \ one corresponding to the initial data and the 
other corresponding to the target data. Points in these manifolds parameter- 
ize corresponding semiclassical states. Then matrix elements of the i?-matrix 
should have the asymptotic 

R«^(u)(a 1 ,<T 2 )= const exp g^lK a ^ ) ^ H {a ll a 2 ){\ + 0(h)) 

h 

where S is the generating function for the mapping p and H is the Hessian of 
S. This asymptotical behavior of the i?-matrix defines the semiclassical limit of 
the partition function and needs further investigation. 



13.3 The large N limit in spin- 1/2 spin chain as a semi- 
classical limit 

13.3.1 

Local spin operators in a spin- 1/2 spin chain of length N are 
S% = 1 <g> . . .< ® • • ■ ® 1, n = l,...,N 

They commute as 



As N — > oo the operators converge to local continuous classical spin 
variables 

77 

s a n -> s%-) 

where S a (x) are local functionals on the classical phase space of the continuum 
spin system. The commutation relations become the relations between Poisson 
brackets: 

{S a (x),S b (y)} = 5(x-y)J2fc b S c (x) 

c 

where S a (x) are continuous local classical spin variables. This combination of 
classical and continuous limit looks more convincing in terms of observables 

1 N 

n=l 
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As N — ¥ oo, such operators becomes S a [f] = J Q f(x)S a (x)dx and the commu- 
tation relation 

[S%[f],S b N [g}] = ^J2 e ^ S »Ud 

C 

becomes 

{S a [flS b [g]} = J2^ b cS c [fd] 

c 

Here we used the asymptotic of commutators in the semiclassical limit: [a, 6] = 
ih{a, b} + . . . and the fact that -h plays the role of the Plank constant. 



13.3.2 

Consider the limit N — > oo of the row-to-row transfer-matrix in the 6- vertex 

2 

model. Assume that at the same time A — > 1 such that A = 1 + + . . . . To 
be specific consider the case when A > 1, so that n = Nr, is finite and real. 
We can write the transfer-matrix of the 6- vertex model as 

t(u) = (smh(u)) N tr(T N (u)) 7 

where T/v(it) is the solution to the difference equation 

T n+ i(u) = Rn }\\ T n (u) 
sinn(it) 

with the initial condition Tq(u) = 1. Here R(u) is the i?-matrix of the 6-vertx 
model. For small r\ we have: 



= 1 + r?(- coth( M )aV^ + " - " ) + 0(r,') 



sinh(u) 2 " sinh(u 

Taking this into account we conclude that as TV — >• 00, T n (u) —> T(u\x) where 
T(u\x) is the solution to 

dT(u\x) A , 1 / \ s^s/ x a + S~(x) + <j~S + (x) 
-UJ. = (-catHuysHx) + ^ ) T i u \ x ) 

with the initial condition T(u\0) — 1. Here a z is the same as in (|32|) . cr ± = 
\(a x ± itrv), and S ± = ^{S 1 ± iS 2 ) 

The row-to-row transfer-matrix of the 6-vertex model have the following 
asymptotic in this limit: 

t(u) -> (sinh(w)) N V(«) (94) 

where r(it) = ir(T(u|l). 
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13.3.3 



Now let us study the evolution of local spin operators in this limit. 

Consider the partition function of the 6- vertex model on a cylinder as a liner 
operator 

^M,JV = K U ) M 

acting in C 2 ®^. Consider the row-to- row transfer-matrix as the evolution op- 

(c) 

erator by one step and Z y M ' N as the evolution operator by M steps. In the 
Heisenberg picture, local spin operators evolve as: 

= t^S^tiu)- 1 (95) 

In the continuum semiclassical limit N —> oo the transfer-matrix becomes a 
functional (|94[) on the phase space of continuous classical spins. The equation 
(|9"5"|) can be written as 

As N — > oo and x — n/N and t — m/N are fixed, this equation become the 
evolution equation for continuum spins: 

d -^± = { H{u) lS «(x,t)} 

where H(u) = \ogt(u). 

The evolution with respect to the partition function on a cylinder of height 

M 

co 7(C) ™ 7(C) 

°n,m+M — Zj M,N°n,m Zj M,N 

becomes the evolution in time T = M /N: S a (x, t) ^ S a (x, t + T). 

If we choose a Lagrangian submanifold in the phase space of continuous spin 
system corresponding to the initial and target data (say, a version of initial and 
target ^-coordinates), the asymptotic of the partition function should be of the 
form 

Z N,m(^^) = conste- NST ^ ^ y/Hess{o x , <r 2 ) (1 + 0(1 /N)) 

where St is the Hamilton- Jacobi action for the Hamiltonian H(u), Hess(a,T) 

(c) 

is the Hessian of St, and in the left side we have the matrix element of Z N M 
between semiclassical states corresponding to o\ and to <t 2 - 

The Bethe equations and the spectrum of the Heisenberg Hamiltonian in 
this limit was studied in [T2l[9l. 



14 The free fermionic point and dimer models 

Decorate the square grid inserting a box with two faces to each vertex, as it is 
shown on Fig. [22l Recall that a dimer configuration on a graph is a perfect 



G9 




CI 



b 



c, 



Figure 21: Projection of dimer configurations on G n m onto the different types 
of vertices at the (n, w)-vertex. 

matching on a set of vertices connected by edges. In other words, it is a collection 
of "occupied edges" (by dimers) such that two occupied edges never meet, and 
any vertex in an endpoint to an occupied edge. 

Dimer configurations on the decorated square grid project to 6- vertex con- 
figurations on a square grid as it is shown on Fig. [5TJ 

It is easy to check that any edge weight system on the decorated lattice 
projects to 6- vertex weights at the free fermionic point, when a 2 + b 2 — c 2 = 
at every vertex. Recall that edge weights is a mapping w : Edges — > R>o- The 
weight of a dimer configuration is 

W{D) = Y[ w(e) 

eG-D 

The statement above means 

E W(D)=Y[w v (S) 
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Figure 22: Reference dimer configuration on G„m- 

where S is a 6- vertex configuration on the square grid, tt is the projection from 
dimer configurations on the decorated square grid to the 6- vertex configurations, 
and the 6- vertex weights w v are given by an explicit formula. The weights w v 
satisfy the free fermionic condition. 

A pair of dimer configurations D, Dq on bipartite graphs define the height 
function. This height function agrees with the height function of the 6-vertex 
model: 

h s (f) = h D , Do (f) 

where S is a 6-vertex state on the square lattice, and h$ is the corresponding 
height function, D 6 7r _1 (S'), Do is shown on Fig. [22j and / is the face of the 
decorated lattice which projects to a face of the square grid. 

References 

[1] Allison, David; Reshetikhin, Nicolai Numerical study of the 6-vertex model 
with domain wall boundary conditions. Ann. Inst. Fourier (Grenoble) 55 
(2005), no. 6, 1847-1869. 

[2] V. I. Arnold. Mathematical Methods of Classical Mechanics, Second Edi- 
tion. Springer, 1989. 

[3] R.J. Baxter, Exactly Solved Models in Statistical Mechanics, Academic 
Press, San Diego, 1982. 

[4] Bogolyubov, N. M.; Izergin, A. G.; Reshetikhin, N. Yu. Finite-size effects 
and critical indices of ID quantum models. JETP Lett. 44 (1986), no. 9, 
521-523 (1987); translated from Pisma Zh. ksper. Teoret. Fiz. 44(1986), 
no. 9, 405-407. 

[5] Bazhanov, V. V.; Reshetikhin, N. Restricted solid-on-solid models con- 
nected with simply laced algebras and conformal field theory. J. Phys. A 
23 (1990), no. 9, 1477-1492. 



71 



D.J. Bukman and J.D. Shore, The Conical Point in the Ferroelectric Six- 
Vertex Model, J. Stat. Phys. 78, 1277-1309, 1995. 

M.Jimbo, T.Miwa, F.Smirnov, Hidden Grassmann Structure in the XXZ 
Model III: Introducing Matsubara direction, arXiv:0811.0439 

V.Chari and A.Pressley, A guide to quantum groups Cambridge University 
Press, 1994. 

T. Bodineau, B. Derrida, V. Lecomte, F. van Wijland Long range 
correlations and phase transition in non-equilibrium diffusive systems, 
larXiv:0807.2394l 

de Vega, H.J., and Destri C, Unified approach to thermodynamic Bethe 
Anstaz and finite size corrections for lattice models and field theories. Nucl. 
Phys. B438, 413-454 (1995). 

V. G. Drinfeld. Quantum groups. In Proc. Intern. Congress of Math. 
(Berkeley 1986), pages 798-820. AMS, 1987. 

Nikolay Gromov, Vladimir Kazakov, Double Scaling and Finite Size Cor- 
rections in sl(2) Spin Chain, Nucl.Phys. B736 (2006) 199-224. 

Skljanin, E. K.; Tahtadjan, L. A.; Faddeev, L. D. Quantum inverse problem 
method. I. (Russian) Teoret. Mat. Fiz. 40 (1979), no. 2, 194-220. 

Tahtadjan, L. A.; Faddeev, L. D. The quantum method for the inverse 
problem and the XYZ Heisenberg model. (Russian) Uspekhi Mat. Nauk 
34 (1979), no. 5(209), 13-63, 256. 

H.Y. Huang, F.Y. Wu, H. Kunz, D. Kim, Interacting Dimers on the Hon- 
eycomb Lattice: an Exact Solution of the Five- Vertex Model, Physica A 
228, 1-32, 1996. 

V. Korepin and P. Zinn- Justin, Thermodynamic Limit of the Six- Vertex 
Model with Domain Wall Boundary Conditions, J. Phys. A, 33, 7053- 
7066, 2000; Inhomogeneous Six- Vertex Model with Domain Wall Boundary 
Conditions and Bethe Ansatz, J. Math. Phys. 43, 3261-3267, 2002. 

M. Jimbo and T. Miwa, Algebraic Analysis of Solvable Lattice Models, 
CBMS Regional Conference Series in Math. 85, 1993. 

Kerov, S. V.; Kirillov, A. N.; Reshetikhin, N. Yu. Combinatorics, the Bethe 
ansatz and representations of the symmetric group. (Russian) Zap. Nauchn. 
Sem. Leningrad. Otdel. Mat. Inst. Steklov. (LOMI) 155 (1986), Differ- 
entsialnaya Geometriya, Gruppy Li i Mekh. VIII, 50-64, 193; translation 
in J. Soviet Math. 41 (1988), no. 2, 916-924. 

[19] V.E. Korepin, N.M. Bogolyubov, and A.G. Izergin, Quantum Inverse Scat- 
tering Method and Correlation Functions, Cambridge University Press, 
1993. 



72 



[20] R. Kenyon, A. Okounkov, and S. Sheffield, Dimers and Amoebae, 
|math-ph/0311005| 

[21] Korogodski, Leonid I.; Soibelman, Yan S. Algebras of functions on quan- 
tum groups. Part I. Mathematical Surveys and Monographs, 56. American 
Mathematical Society, Providence, RI, 1998. 

[22] Kulish, P. P.; Reshetikhin, N. Yu.; Sklyanin, E. K. Yang-Baxter equations 
and representation theory. I. Lett. Math. Phys. 5 (1981), no. 5, 393-403. 

[23] E. Licb, Phys. Rev. 162, 162, 1967; E.H. Lieb, Phys. Rev. Lett. 18, 1046, 
1967; 19, 108, 1967. 

[24] E.H. Lieb and F.Y. Wu, Two Dimensional Ferroelectric Models, in Phase 
Transitions and Critical Phenomena, Vol. 1, ed. by C.Domb and M.S. 
Green, 321, Academic Press, London, 1972. 

[25] B. Nienhuis, H.J. Hilhorst, and H.W. Bloete, Triangular SOS modles and 
cubic-crystal shapes, J. Phys. A: Math. Gen., 17, (1984), 3559-3581. 

[26] J.D. Noh and D. Kim, Finite-Size Scaling and the Toroidal Partition Func- 
tion of the Critical Asymmetric Six- Vertex Model, |cond-mat / 95 1 1001} 

[27] I. M. Nolden, The Asymmetric Six- Vertex Model, J. Statist. Phys. 67, 155, 
1992; Ph.D. thesis, University of Utrecht, 1990. 

[28] K. Palamarchuk, Ph.D. thesis, Berkeley, 2007. 

[29] K. Palamarchuk, and N. Reshetikhin, The 6- 

vertex model with fixed boundary conditions, 
http://pos.sissa.it/archive/conferences/038/012/Solvay-012.pdf. 

[30] Reshetikhin, N. Quasitriangularity of quantum groups at roots of 1. Comm. 
Math. Phys. 170 (1995), no. 1, 79-99. 

[31] Reshetikhin, N. Yu.; Smirnov, F. A. Quantum Floquet functions. (Russian) 
Questions in quantum field theory and statistical physics, 4. Zap. Nauchn. 
Sem. Leningrad. Otdel. Mat. Inst. Steklov. (LOMI) 131 (1983), 128-141. 

[32] Reshetikhin, Nicolai; Yakimov, Milen Deformation quantization of La- 
grangian fiber bundles. Confrence Mosh Flato 1999, Vol. II (Dijon), 263- 
287, Math. Phys. Stud., 22, Kluwer Acad. PubL, Dordrecht, 2000. 

[33] A.Reiman, M.Semenov-Tian-Shanskii, Group theoretical methods in In- 
tegrable systems, in: Integrable systems VII (V.Arnold and S. Novikov, 
eds.) Encyclopaedia of Math Sciences v. 16 Springer 1994 

[34] Skljanin, E. K. The method of the inverse scattering problem and the quan- 
tum nonlinear Schrdinger equation. (Russian) Dokl. Akad. Nauk SSSR 244 
(1979), no. 6, 1337-1341. 



73 



[35] Sklyanin, E. K. Classical limits of SU(2)-invariant solutions of the Yang- 
Baxter equation. (Russian) Translated in J. Soviet Math. 40 (1988), no. 1, 
93-107. Differential geometry, Lie groups and mechanics. VII. Zap. Nauchn. 
Sem. Leningrad. Otdel. Mat. Inst. Steklov. (LOMI) 146 (1985), 119-136, 
203, 206-207. 

[36] Smirnov, F. A. Quasi-classical study of form factors in finite volume. L. D. 
Faddeev's Seminar on Mathematical Physics, 283-307, Amer. Math. Soc. 
Transl. Scr. 2, 201, Amer. Math. Soc, Providence, RI, 2000. 

[37] B. Sutherland, C.N. Yang, and CP. Yang, Exact Solution of a Model of 
Two-Dimcnsional Ferroelectrics in an Arbitrary External Electric Field, 
Phys. Rev. Letters 19, 588, 1967. 

[38] Yakimov, Milen Symplectic leaves of complex reductive Poisson-Lie groups. 
Duke Math. J. 112 (2002), no. 3, 453-509. 

[39] C.N.Yang and C.P.Yang Thermodynamics of a one- dimensional system of 
bosons with repulsive delta-function interactions, J. Math. Phys., 10 (1967) 
1115-1122 

[40] Al. Zamolodchikov, Thermodynamic Bethe anstz in relativistic models: 
Scaling 3-state-Potts and Lee- Yang models. Nucl. Phys, B342, 695-720 
(1990). 

A Symplectic and Poisson manifolds 
A.l 

Recall that an even dimensional manifold equipped with a closed non-degenerate 
2-form is called symplectic. 

Let (M, to) be a 2n-dimensional symplectic manifold. In local coordinates: 

In In 2n q 

lu = ^2 Uijdx 1 A dx\ det(u>) ^ 0, duj = ^ -^^dx 1 A dx' J = 0, 

ij=l fe=l ij = l 

The last identity is equivalent to the Jacobi identity for the bracket 
A.2 

A smooth manifold M with bi- vector field p (a section of the bundle A 2 TM) 
such that the bracket between two smooth functions 

{f,9}=p(df Adg) 

satisfies the Jacobi identity is called a Poisson manifold. I local coordinates 
x 1 , . . . , x n , p(x) = £ i, j = 1 V W^A^. 
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A.3 



A Poisson tensor on a smooth manifold M defines a subspace mapping p : 
T*M -> TM. The image is a system of subspaces p(T*M) C TM which is a 
distribution on M. Leaves of this distribution are spanned by curves which are 
flow lines of piece- wise Hamiltonian vector fields. They are smooth submani- 
folds. Symplectic leaves of the Poisson manifold M are leaves of this distribution. 

Well know examples of symplectic leaves are co-adjoint orbits in the dual 
space to a Lie algebra. 

B Classical integrable systems and their quan- 
tization 

B.l Integrable systems in Hamiltonian mechanics 

The notion of integrability is most natural in the Hamiltonian formulation of 
classical mechanics. For details see [2] [33]. 

In Hamiltonian formalism of classical mechanics the dynamics is taking place 
on the phase space and equations of motion are of the first order. When it is a 
system of particles moving on a manifold M, local coordinates are positions and 
momenta of particles. Globally, the phase space in this case is the cotangent 
bundle to the manifold M, see for example [2]. 

Hamiltonian formulation of spinning tops or other systems with more com- 
plicated constraints involves more complicated phase spaces. In all cases a phase 
space has a structure of a structure of a symplectic manifold, that is it comes 
together with a non-degenerate closed 2-form on it. 

Any symplectic manifold (and so any phase space of a Hamiltonian system 
admits local coordinates (Darboux coordinates) in which the 2-form has the 
form 

n 

to = dpi A dq l 
»=l 

These coordinates can be interpreted as momenta and positions, though in case 
of spinning tops this interpretation does not have a lot of physical meaning but 
Darboux coordinates is a convenient mathematical tool. 

The dynamics in a Hamiltonian system is determined by the energy function 
H. The trajectories of such system in local Darboux coordinates are solutions 
to differential equations: 

dq % _ dH d Pi _ dH 
dt dpi ' dt dq l 

Geometrically, trajectories are flow lines of the Hamiltonian vector field vh = 
u~ l {dH) G r(A 2 TM) where a;" 1 : TM -t T*M is the bundle isomorphism 
induced by the symplectic form lu. 

Let M be a 2iV-dimensional symplectic manifold. 
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Definition 1. An integrable system on M is a collection of N independent 
functions on M which commute with respect to the Poisson bracket. 

Recall that 



• The "level surfaces' 



M(ci,...,c„) = {xe M,Ii(x) =a} (96) 
are invariant with respect to the flow of any Hamiltonian H = F(Ii, . . . ,In)- 

• For every such Hamiltonian and every level surface M(ci, . . . , c„) there 
exists an affine coordinate system (p±, . . . ,p n ) in which the Hamiltonian 
flow generated by H is linear: pi =constant. 

• The coordinate system (jpi, . . . ,p n ) on M(cx, . . . , c n ) can be completed to 
a canonical coordinate system (pi, qj) in every sufficiently small neighbor- 
hood of M(cx, . . . , c„). These coordinates are called action- angle variables 
and in some cases these coordinates are global. 

C Poisson Lie groups 

C.l Poisson Lie groups 

A Lie group G is called a Poisson Lie group if 

• G has a Poisson structure, i. e. it is given together with the Poisson 
tensor p € T(A 2 TG), in local coordinates p(x) = p ij (x) A -£j <G A 2 T X G. 
This tensor the Poisson bracket on smooth functions on G (Lie bracket 
satisfying the Leibnitz rule with respect to the point-wise multiplication) 

{/,.}(*) = £p«(*)g||(z) 



• The Poisson structure is compatible with the group multiplication. This 
means that the multiplication mapping G x G — > G bring the Poisson 
tensor on G x G to the Poisson tensor on G, or: 

T.p l3 ^y)^^ z )\^y - Ep ij ( x )^^(xy)+J2p ij (y)^^(*y) 

ij ij ij 

for any pair of functions /, g. 
For more details on Poisson Lie groups and for numerous examples see 
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C.2 Factorizable Poisson Lie groups and classical r-matrices 



The class of Poisson Lie groups relevant to integrable systems has so-called r- 
matrix Poisson brackets. Let g be a Lie algebra corresponding to the Lie group 
G. A classical r-matrix for G is an element r £ g (g> g satisfying the bilinear 
identity 

[ri2,r 13 ] + [ri2, r 23 ] + [r 13 , r 23 ] = (97) 

The Poisson tensor 

p(x) =r~ Ad x (r) £ g A g 

defines a Poisson Lie structure on G if r satisfies (|97fl and r + a(r) is an invariant 
tensor. 

Two remarks: wc identified T x GwithT e G = g using left translations on G; 
Ad x is the diagonal adjoint action of G on g A g, u A v — > xux~ 1 ® xvx~ x when 
G is a matrix Lie algebra. 

Let 7r y : G -> G£(V) and tt^ : G -> GL(M / ) be two (finite-dimensional) 
representations of G and let 7r^,7r^ be matrix elements of G in a linear basis: 
7Ty : g € G — >• TrYj(g) £ C. The r-matrix Poisson brackets between two such 
functions are: 

Taking trace in this formula we see that characters of finite dimensional represen- 
tations form a Poisson commutative algebra on G. Restricting this subalgebra 
to a symplectic leaf of G we will obtain an integrable system if the number of 
independent functions among \\ after this restriction is equal to half of the 
dimension of the symplectic leaf. 

the principal advantage of this approach is that it gives an algebraic way to 
construct classical r-matrices through the double construction of Lie bi-algebras. 

In then next section we will see how r-matrices with spectral parameter 
appear naturally from Lie bi-algebras on loop algebras. 

C.3 Basic example G = LSL 2 

Our basic example is an infinite dimensional Poisson Lie group LGLi of map- 
pings S 1 = {z £ C\\z\ = 1} — > GL2 which are holomorphic inside the unit 
disc. 

The Lie algebra Lgl 2 (we consider maps which are Laurent polynomials in 
t) of the Lie group LGL 2 has a linear basis -[n] with ey [n](z) = eyt - ™ . It 
also has an invariant scalar product (x[n},y[m\) = 8 n .m- The following element 
of the completed tensor product Lgl 2 ® Lgl 2 

r= ^ eij [0] ® eji [0] + ^ e ij[ n ] ® eji[-n] (98) 

i>j;i,i=l,2 i j' = l,2;neZ,ra>l 

satisfies the classical Yang-Baxter equation (|97"1) . Wc will not explain here why, 
but this follows from the Drinfeld's double construction for Lie bialgebras. 
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Let tt v : GL-2 — > GL(V) be a representation of GL2 and a 6 C v be a non- 
zero complex number. The evaluation representation ny^ a of Lg^ acts in V 
as 

7r Vi A4n])=7r v (x)a- 2n - 2 

where a is a non-zero complex number. 

Evaluating the element (|98|) in iry.a ® 7i"v,6 where V is the two dimensional 
representation we get 

nv,a ® 7ry,6(r) = r(a/i>) + f(a/b)l 

where /(z) is a scalar function and 

z + z -1 2 

r(z) = — T a 7 -®a z A -(cr+ ® <j~ + a~ ® cr+) (99) 

z — z 1 z — z 1 

Here <r z = en — e22, cr + = £12, and cr _ = e2i are Pauli matrices. 
The Yang-Baxter equation for r implies 

[r 12(2), r 13 (zw)} + [r 12 {z),r 23 (w)} + [ri 3 (zw), r 2 3(w)] = 

The Poisson bracket between coordinate functions g%j{z) are 

{ gi (z),g 2 (w)} = [r(z/w), gi (z)g 2 (w)} (100) 

One of the important properties of such Poisson brackets is that charac- 
ters of finite dimensional representations define families of Poisson commuting 
functions on LG2- 

t v (z) = tr v (ir v (g(x)), {t v (z) , t w (w)} = 

The coefficients of these functions will produce Poison commuting integrals of 
motions for integrable systems, when restricted to symplectic leaves of LGL2- 

C.4 Symplectic leaves 

It is a well know, classical fact, which can be traced back to works of Lie, that 
symplectic leaves of the Poisson manifold which is the dual space to a Lie algebra 
are co-adjoint orbits. Similarly, symplectic leaves of a Poisson Lie group G are 
orbits of the dressing action of the dual Poison Lie group on G. 

The structure of symplectic leaves of Poisson Lie groups is well known for 
finite dimensional simple Lie algebras, see [H] 38 . For the construction of in- 
tegrable spin chains related to SL2 we will need only some special symplectic 
leaves of LGL2- These symplectic leaves are symplectic leaves of finite dimen- 
sional Poisson submanifolds of polynomial maps of given degree. 
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D Quantization 



D.l Quantization 

This section is a brief outline of the quantization of Hamiltonian systems. There 
is also a very important point of view of a path integral quantization. We will 
not discuss it here. 

D.l.l Quantized algebra of observables 

Let M be a symplcctic manifold, which is the phase space of our mechanical 
system. We want to describe possible quantum mechanical systems which re- 
produce our system in the classical limit. This procedure is called quantization. 

Let A be a Poisson algebra over C i.e. it is a complex vector space with 
a commutative multiplication ab and with the Lie bracket {a, b} such that 
{a, be} — b{a, c} + {a, b}c. Let X C C is a neighborhood of e C. 

Definition 2. A deformation quantization of A is a family of associative al- 
gebras Ah, parameterized by h e X together with two families of linear maps 
4>h : Ah — > A and iph '■ A — > Ah such that 

• lim^o (ph°^h^ id,A, 

• lim/,->o 4>h(iph(a)iph(b)) = ab, 

• lim^o ih(M^M^M2^thM = { a , & } ; 

Denote by C(M) the classical algebra of observables, i.e. the algebra of real 
valued functions on the phase space. If M is T*N, this will be the algebra of 
functions on M which are polynomial in the cotangent direction and smooth on 
N. If M is an affine algebraic, manifold, C(M) will be algebra of polynomial 
functions. If M is a smooth manifold it is C°°M, etc.. The space C(M) has 
a natural structure of a Poisson algebra with the point-wise multiplication and 
the Poisson bracket determined by the symplectic form on M. 

Denote its complcxification by C(M)c- The real subalgebra C(M) C C(M)c 
is the set of fixed points of complex conjugation, i.e. real valued functions. 
Denote complex conjugation by a, i. e. a(f )(x) = f{x). 

Deformation quantizations which are relevant to quantum mechanics are real 
deformation quantizations: 

Definition 3. A real deformation quantization ofC(M) is a pair (Ah,o~h,h G 
X C R) where Ah is a deformation quantization (as above) of C(M)c, h is a 
real deformation parameter, and ah is a C- anti-linear anti-involution of Ah, i.e. 

o-l = l, o- h (ab) = a h (b)a h (a), a h (sa) = sa h {a), 

The the subspace C{M)h C Ah of fixed points of ah is called the space of quantum 
observables. Linear maps 4>h and iph from the definition f Ah should satisfy extra 
condition 

lim 4>h o ah o iph = & 

h->0 
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The space of fixed points of ah is called the space of quantum observables. 
Notice that the multiplication does not preserve this space. However, it is closed 
with respect to operations AB + BA and i(AB — BA). Such structures are called 
Jordan algebras. 

Of course these definition still need a clarification. The vector space C (M) 
is infinite dimensional and we have to specify in which topology our linear 
isomorphism is continuous and in which sense we should take limits. The way 
how to handle this is either dictated by the nature of a specific problem, so 
we will discuss it later in relation to integrable spin chains which is the main 
subject here. 

D.2 Examples of family deformations 
D.2.1 

Take M = R 2 , A = Pol c (M, 2 ) = C\p,q] with the standard symplectic form 
dp A dq giving the bracket {p, q} = 1 (this determines the bracket). We have a 
natural monomial basis p n q m on A. Define 

A h =< p, q\pq -qp = h> 

It is clear that this is a family of associative algebras. To identify this family 
with a deformation quantization of A we should find <f>h and iph- For this choose 
monomial bases p n q m in Ah and A. Define linear maps <f>h and iph as linear 
isomorphisms Ah — A identifying the monomial bases. It is clear that this 
choice makes Ah into a defromation quantization of A. 

D.2. 2 

Let g be a Lie algebra, and consider Pol(g*) = C[g]. If {e^} is a basis for g, 
then we can think of the e, as coordinate functions xi on g*. A theorem of 
Kostant, {f,g}(x) =< x,[df(x),dg(x)] >, {xi,xj} = Y.k c % x k- Pol(g*) then 
gets a Poisson bracket. 

We can get a deformation quantization 

Ah ^ X\ , . . . , X n \ XiXj XjX{ H ^ ^ C^jXfc ^> 

k 

Note that Ah = Ug for any h ^ (you just have to rescale the x's by h). On 
the other hand, Choosing the monomial basis 1 and 

the PBW basis in Ah- Identifying them, we get Ah = C[x\, . . . ,x n ], which is 
how the PBW theorem is usually formulated. 

Ug = Pol(g*) ^ S(g) 

We get a linear isomorphism <f>h : Ah = A. It is easy to check that this is a 
deformation quantization. 
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D.3 Quantization of integrable systems 

As we have seen above the quantization of a classical Hamiltonian system on M 
with the Hamiltonian H £ C(M)consists of the following: 

• A family of associative algebras Ah quantizing the algebra of function on 
M. 

• The choice of quantum Hamiltonian Hh £ Ah for each /i,such that <f>h {Hh) ~> 
H as h — > 0. 

The quantization is integrable if for each h there is a maximal commutative 
subalgebra of Ch(M) quantizing the subalgebra of classical integrals in C(M), 
which contains Hh- 

To be more precise, let I(M) £ C{M) be the subalgebra generated by Pois- 
son commuting integrals in the classical algebra of observables, i.e. {H, F} = 
for each F £ I(M) and {G, F} = for each F,G £ I(M). Its integrable quan- 
tization consists of a commutative subalgebra Ih £ Ah such that FH — HF for 
each F £ Ih, such that lim^o <fih(Ih) — I- For the precise definition in case of 
formal deformation quantization see [32] . 

When quantized algebra of observables is represented in a Hilbert space, the 
important problem is the computation of the spectrum of commuting Hamilto- 
nians. 
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